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Abstract 

^ ■ 

. We consider the transverse-momentum (gy) distribution of generic high-mass systems 

• ^ ! (lepton pairs, vector bosons, Higgs particles, ....) produced in hadron collisions. 

^ I At small qx, we concentrate on the all-order resummation of the logarithmically- 

' enhanced contributions in QCD perturbation theory. We elaborate on the 6-space 

resummation formalism and introduce some novel features: the large logarithmic 
contributions are systematically exponentiated in a process-independent form and, 
after integration over qt, they are constrained by perturbative unitarity to give a 
vanishing contribution to the total cross section. At intermediate and large qt, 
resummation is consistently combined with fixed-order perturbative results, to obtain 
predictions with uniform theoretical accuracy over the entire range of transverse 
momenta. The formalism is applied to Standard Model Higgs boson production at 
LHC energies. We combine the most advanced perturbative information available 
at present for this process: resummation up to next-to-next-to-leading logarithmic 
accuracy and fixed-order perturbation theory up to next-to-leading order. The results 
show a high stability with respect to perturbative QCD uncertainties. 
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1 Introduction 



This paper is devoted to study the transverse-momentum [qt) spectrum of high-mass systems 
produced by hard-scattering of partons in hadron-hadron colhsions. In Ref. P we presented some 
quantitative results on the qt spectrum of the Standard Model (SM) Higgs boson, produced via 
the gluon fusion mechanism, at LHC energies. The formalism used in Ref. P is quite general 
and applies to the transverse-momentum distribution of generic high-mass systems (lepton pairs, 
vector bosons, Higgs particles, ..) produced in hadron collisions. The purpose of the present paper 
is twofold. Owing to its general applicability, we find it useful to first describe and discuss the 
formalism with quite some details. We then perform a more systematic phenomenological analysis 
of the qt distribution of the Higgs boson at the LHC. 

In this introductory section, rather than illustrating the resummation formalism in general 
terms, we mainly consider the explicit case of the qt spectrum of the Higgs boson. This also 
serves for underlying some general features of the formalism in concrete, rather than abstract, 
terms. 

Within the SM of electroweak interactions, the Higgs boson [2j is responsible for the mechanism 
of the electroweak symmetry breaking, but this particle has so far eluded experimental discovery. 
Direct searches at LEP have established a lower bound of 114.4 GeV [3j on the mass Mh of the SM 
Higgs boson, whereas SM fits of electroweak precision data lead to the upper limit Mjj < 260 GeV 
at 95% CL The next search for Higgs boson(s) will be carried out at hadron colliders, namely, 
the Fermilab Tevatron 0^1 and the CERN LHC [71IH]. 

The main production mechanism of the SM Higgs boson H at hadron colliders is the gluon 
fusion process gg —>■ H, through a heavy-quark (mainly, top-quark) loop. When combined with the 
decay channels H — > 77, H — > WW and H ZZ, this production mechanism is one of the most 
important for Higgs boson searches and studies over the entire mass range, 100 GeV < ^ 1 TeV, 
to be investigated at the LHC [7j. To fully exploit the physics potential of the gluon fusion process, 
it is relevant to provide reliable theoretical predictions for the corresponding total cross section 
and for the associated distributions, such as, for instance, the Higgs qx distribution. The dominant 
source of theoretical uncertainties on these quantities is the effect of QCD radiative corrections, 
which, therefore, have to be carefully investigated. 

The total cross section for Higgs boson production by gluon fusion has been computed in QCD 
perturbation theory at the leading order (LO), C(ag), at the next-to-leading order (NLO) P fTU] 
and at the next-to- next-to-leading order (NNLO) [TTHT^ in the QCD coupling ag. The NNLO 
computation of the rapidity distribution of the Higgs boson has recently been completed ^H]- A 
key point of this theoretical activity is that the origin of the dominant perturbative contributions 
to the total cross section has been identified and understood: the bulk of the radiative corrections 
is due to virtual and soft-gluon terms [12]. This point has a twofold relevance. On one side, it 
explains the observation of the validity of the large-M^ approximation (Mf being the mass 
of the top quark) in the calculation at the NLO, and, therefore, it justifies the use of the same 
approximation at and beyond the NNLO. On the other side, it allows to estimate higher-order 
QCD contributions by supplementing the NNLO calculation with an all-order resummation of 
the logarithmically-enhanced terms due to multiple soft-gluon emission ^7]. Having these terms 
under control allows us to reliably predict the value of the cross section and, more importantly. 
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to reduce the associated perturbative uncertainty at the level of about ±10% [T7j . 

When studying the qx distribution of the Higgs boson in QCD perturbation theory, it is 
convenient to start by considering separately the large-g^ and small-gr regions. 

The large-gT region is identified by the condition qx ~ M^. In this region, the perturbative 
series is controlled by a small expansion parameter, as(M|f), and calculations based on the trun- 
cation of the series at a fixed order in as are theoretically justified. SM Iliggs boson production at 
large via gluon fusion has to be accompanied by the radiation of at least one recoiling parton, 
so the LO term for this observable is of 0{a\). The LO calculation was reported in Ref. |18j : 
it shows that the large- Mt approximation works well as long as Mn^^Mt and qr^Mt. Similar 
results on the validity of the large- approximation were obtained in the case of the associated 
production of a Higgs boson plus 2 jets (2 recoiling partons at large transverse momenta) 
In the framework of the large-M^ approximation, the NLO QCD corrections to the transverse- 
momentum distribution of the SM Higgs boson were computed in Refs. pU^I23j . Corrections to 
the large- approximation are considered in Ref. [21] . The numerical programs of Refs. [201 122] 
can also be used to evaluate arbitrary infrared- and coUinear-safe observables up to NLO in the 
large-gr region and, in the case of Ref. |22]; up to NNLO when q^ = 0. 

In the small-gr region {q^ ^ M^), where the bulk of events is produced, the convergence of 
the fixed-order expansion is spoiled, since the coefficients of the perturbative series in as{M]j) 
are enhanced by powers of large logarithmic terms, In^ {Mfj / q^) . To obtain reliable perturbative 
predictions, these terms have to be resummed to all orders in as. The method to systematically 
perform all-order resummation of classes of logarithmically-enhanced terms at small gr is known 
[2SHSS]- In the case of the SM Higgs boson, resummation has been explicitly worked out at 
leading logarithmic (LL), next-to- leading logarithmic (NLL) IHS] and next-to- next-to-leading 
logarithmic (NNLL) ^ level. 

The fixed-order and resummed approaches at small and large values of gr can then be matched 
at intermediate values of qr, to obtain QCD predictions for the entire range of transverse momenta. 
Phenomenological studies of the SM Higgs boson q^ distribution have been performed in Refs. [23 

I^HH^ in E3ti46j . by combining resummed and fixed-order perturbation theory at different levels 
of theoretical accuracy. A comparison of theoretical calculations ^ HHl S21 and of results from 
parton shower Monte Carlo generators [TfHSU] is presented in Ref. |51j . 

In the present paper we compute the Higgs boson distribution at the LHC by combining the 
most advanced perturbative information that is available at present: NNLL resummation at small 
g-r and NLO perturbation theory at large g^. The first results of our calculation were presented in 
Refs. ^ 122] • Here we perform a more complete phenomenological study and present a discussion 
of theoretical uncertainties. 

The formalism used to obtain these results was briefiy described in Refs. [321 d] and is illustrated 
in detail in the present paper. Three distinctive features are anticipated here. The resummation 
is performed at the level of the partonic cross section; this implies that the parton distributions 
are evaluated at the factorization scale fip, which has to be chosen of the order of the hard scale 
M. The resummed terms are embodied in a form factor that is universal: it depends only on 
the fiavour of the partons that initiate the hard-scattering subprocess at the Born level (e.g. gg 
annihilation in the case of Drell-Yan lepton pair production, and gg fusion in the case of Higgs 
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boson production). A constraint of perturbative unitarity is imposed on the resummed terms, to 
the purpose of reducing the effect of unjustified higher-order contributions at large values of qx 
and, especially, at intermediate values of qx- The constraint implies that the total cross section 
at the nominal fixed-order accuracy (NLO or NNLO) is recovered upon integration over qx of the 
transverse- momentum spectrum (at NLL+LO or NNLL+NLO). 

The paper is organized as follows. In Sect.Elthe resummation formalism is discussed in detail. 
After illustrating the general aspects of our approach in Sect. 12.11 we discuss the structure of the 
resummed cross section in Sect. 12.21 The relation to the standard 6-space resummation is given 
in Sect. 12.31 Sect. 12.41 is devoted to the finite component of the cross section. In Sect. |3]we apply 
the resummation formalism to the production of the SM Higgs boson at the LHC In Sect. |3]we 
draw our conclusions. In Appendix ^ we discuss the details of the exponentiation in the general 
multiflavour case. In Appendix^ we illustrate the calculation of the Bessel integrals required in 
the computation of the perturbative expansion of the resummed cross section. 



2 Transverse- momentum resummation 

The formalism [H inB] that we use to compute the qx distribution of the Higgs boson applies to 
more general hard-scattering processes. Therefore, we describe it in general terms. 

2.1 The resummation formalism: from small to large values of qt 

We consider the inclusive hard-scattering process 

h^{pi) + h2{p2)^F{M,qT)+X , (1) 

where the collision of the two hadrons hi and /i2 with momenta pi and p2 produces the triggered 
final-state system F, accompanied by an arbitrary and undetected final state X. We denote by 
\fs the centre-of-mass energy of the colliding hadrons (s = {pi +P2Y — 2piP2)- The observed final 
state F is a generic system of non-QCD partons such as one or more vector bosons (7*, W, Z, . . .), 
Higgs particles, Drell-Yan (DY) lepton pairs and so forth. We do not consider the production of 
strongly interacting particles (hadrons, jets, heavy quarks, ...), since in this case the resummation 
formalism of small-gy logarithms has not yet been fully developed. 

Throughout the paper we limit ourselves to considering the case in which only the total in- 
variant mass M and transverse momentum q-r of the system F are measured. According to the 
QCD factorization theorem (see Ref. [^2] and references therein), the corresponding transverse- 
momentum differential cross section^ dap/dq^ can be written as 




(2) 

^To be precise, when the system F is not a single on-shell particle of mass M, what we denote by dap/dq^ is 
actually the differential cross section M'^dap /dM'^dq^. 
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where fa/hi^yf^p) (f* = Qf, Qf:9) cire the parton densities of the colhding hadrons at the factoriza- 
tion scale lip, dapab/dq^ are the partonic cross sections, s = X1X2S is the partonic centre-of-mass 
energy, and iJ,ji is the renormahzation scale. Throughout the paper we use parton densities as 
defined in the MS factorization scheme, and as(g^) is the QCD running coupling in the MS 
renormahzation scheme. 

The partonic cross section is computable in QCD perturbation theory as a power series ex- 
pansion in as. We assume that at the parton level the system F is produced with vanishing qp 
(i.e. with no accompanying final-state radiation) in the lowest-order approximation, so that the 
corresponding cross section is da^l^/dq^ oc S{q^). Since F is colourless, the lowest-order partonic 
subprocess, c + c — >• F, is either qq annihilation {c = q), as in the case of j*,W and Z production, 
or gg fusion {c = g), as in the case of the production of the SM Higgs boson H. 

As recalled in Sect. ^ higher-order perturbative contributions to the partonic cross section 
dcTpab/dq^ contain logarithmic terms of the type ln™'(M^/g|.) that become large in the small-g^ 
region {qp -C M). Therefore, we introduce the following decomposition of the partonic cross 
section in Eq. (j21): 

7^ 7^ (res.) 7^ (fin.) 

dq^ dq^ dq^ 

The distinction between the two terms on the right-hand side is purely theoretical. The first 
term, da^p^^ , on the right-hand side contains all the logarithmically-enhanced contributions, 
(og/gl^) ln'^(M^/g|,), at small qp, and has to be evaluated by resumming them to all orders 
in as. The second term, da^p^^ , is free of such contributions, and can be computed by fixed-order 
truncation of the perturbative series. More precisely, we define the 'finite' component da^p^'i^ in 
such a way that we have* 

liin / dq^ =0, (4) 



^0 



dq, 



f.o. 



where the right-hand side vanishes order- by- order in perturbation theory. In particular, this 

.(fin.) 
'Fab 



implies that any perturbative contributions proportional to S{q^) have been removed from daf^^ '' 
and included in dap^^^\ 

The 'resummed' component d&p^^f^^ of the partonic cross section cannot, of course, be evaluated 
by computing all the logarithmic contributions in the perturbative series. However, as discussed 
in Sect. 12.21 these contributions can systematically be organized in classes of LL, NLL, . . . terms 
and, then, this logarithmic expansion can be truncated at a given logarithmic accuracy. 

In summary, the qp distribution in Eq. Q is evaluated, in practice, by replacing the partonic 
cross section on the right-hand side as follows 



dq? 



T 



d^Fab 
,2 



dqp J i.a. 



+ 



J- (fin.) 
d(^Fab 

dqp J f.o. 



(5) 



The first and second terms on the right-hand side denote the truncation of the resummed and 
finite components at a given logarithmic accuracy and at a given fixed order, respectively. The 
resummed component gives the dominant contribution in the small-gr region, while the finite 



*The notation [X]^^ means that the quantity X is computed by truncating its perturbative expansion at a 
given fixed order in ag- 
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component dominates at large values of qx- The two components have to be consistently matched 
at intermediate values of g^? so as to obtain a theoretical prediction with uniform formal accuracy 
over the entire range of qr-, from -C M up to qx ~ M. To this aim, we compute ['^'^"^b'^ '']^^ 
starting from \daab\ ^ ^ , the usual perturbative series for the partonic cross section truncated at 
a given fixed order in as, and subtracting from it the perturbative truncation of the resummed 
component at the same fixed order in as'- 



da 



(fin.) 
Fab 



dq. 



f.o. 



da 



Fab 



- dq, 



f.o. 



da 



(res.) 
Fab 



- dq, 



J f.o. 



(6) 



Moreover, we impose the condition: 



7 " (res.) 






7 ^ (res.) 
\d<^Fab^ 




- dq^ . 


La. 


f.o. 


. dq^ . 


f.o. 



(7) 



This matching procedure guarantees that the replacement in Eq. (0) retains the full information 
of the perturbative calculation up to the specified fixed order plus resummation of logarithmically- 
enhanced contributions from higher orders. Equations © and ((7j) indeed imply that the matching 
is perturbatively exact, in the sense that the fixed-order truncation of the right-hand side of 
Eq. (0) exactly reproduces the customary fixed-order truncation of the partonic cross section in 
Eq. (0). The (small-gr) resummed and (large-gr) fixed-order approaches are thus consistently 
combined without double-counting (or neglecting) of perturbative contributions and by avoiding 
the introduction of ad-hoc boundaries (such as, for instance, the choice of some intermediate 
value of qT as 'switching' point between the resummed and fixed-order calculations) between the 
large-gr and small-g^ regions. 

The resummed contributions that are present in the term ['^<5"i7'ab'']i a ® necessary 

and fully justified at small qx- Nonetheless they can lead to sizeable higher-order perturbative 
effects also at large qx, where the small-gr logarithmic approximation is not valid. To reduce the 
impact of these unjustified higher-order terms, we require that they give no contributions to the 
most basic quantity, namely the total cross section, that is not affected by small-gr logarithmic 
terms. We thus impose that the integral over q^ of Eq. © exactly reproduces the fixed-order 
calculation of the total cross section. Since da^p^^ is evaluated in fixed-order perturbation the- 
ory, the perturbative constraint on the total cross section is achieved by imposing the following 
condition: 



dqj 



da 



(res.)_, 
Fab 



dq. 



l.a. 



dq. 



T 



da 



(res.) 
Fab 



dq; 



f.o. 



(8) 



Equation (jSJ can be regarded, in some sense, as a unitarity constraint. As a matter of fact, the 
logarithmic contributions that are resummed in da^p^j^ are, precisely speaking, plus distributions 
of the type [{o-s/It) 1'^'^(^^/5't)] + - Therefore, it is quite natural to require that these resummed 
terms give a vanishing contribution to the total cross section. Note that the bulk of the qr dis- 
tribution is in the region qx ^Mh- Since resummed and fixed-order perturbation theory controls 
the small-gr and large-g^ regions respectively, the total cross section constraint mainly acts on the 
size of the higher-order contributions introduced in the intermediate-g^ region by the matching 
procedure. 

Another distinctive feature of the formalism illustrated so far is that we implement perturbative 
QCD resummation at the level of the partonic cross section. In the factorization formula (j2I), the 
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parton densities are thus evaluated at the factorization scale /j^p, as in the customary perturbative 
calculations at large qt- Although we are dealing with a process characterized by two distinct hard 
scales, qT and M, the dominant effects from the scale region <^ M are explicitly taken into 
account through all-order resummation. Therefore, the central value oi hf and has to be set 
equal to Mh-, the 'remaining' typical hard scale of the process. Then the theoretical accuracy of 
the resummed calculation can be investigated as in customary fixed-order calculations, by varying 
jip and hr around this central value. 

At small values of qx, the perturbative QCD approach has to be supplemented with non- 
perturbative contributions, since they become relevant as qp decreases. A discussion on non- 
perturbative effects on the g-r distribution of the SM Higgs boson is presented in Sect. 13.11 

The resummation and matching formalism, which we have so far illustrated in quite general 
terms, is set up to deal with the transverse-momentum region where qp ^ M. Resummation of 
small-gr logarithms cannot lead to any theoretical improvements in the large-g^ region, where 
those logarithms are not the dominant contributions. When qp ^ M, the use of the resummation 
formalism is no longer justified (recommended), and we have to use the customary fixed-order 
perturbative expansion. 



2.2 The resummed component 

The method to systematically resum the logarithmically-enhanced contributions at small qx was 
set up HHHSOl shortly after the first resummed calculation of the DY qp spectrum to double 
logarithmic accuracy [23- The resummation procedure has to be carried out in the impact- 
parameter space, to correctly take into account the kinematics constraint of transverse-momentum 
conservation. The resummed component of the transverse-momentum cross section in Eq. (jH)) is 
then obtained by performing the inverse Fourier (Bessel) transformation with respect to the impact 
parameter h. We write^ 

7^ (res.) , j-2 r j2i_, 

^^(gT,M,5;as(/x?j),/iL/^F) = — j — e^^"'- >Vf,(6, M, S; as(/i?j), /i^ /^f) (9) 

b 

= — / db-MbqT)W^,{b,M,s;as{fil),fil,fil),m 
s Jo ^ 

where J^ip^^ is the Oth-order Bessel function. 

The perturbative and process-dependent factor embodies the all-order dependence on the 
large logarithms In M^b^ at large b, which correspond to the gr-space terms In M^/g|. that are 
logarithmically enhanced at small gr (the limit -C M corresponds to Mb 3> 1, since b is the 
variable conjugate to g^). Resummation of these large logarithms is better expressed by defining 
the iV-moments^ Wn of W with respect to z = M'^/s at fixed M: 

W^,^j,{b,M;as{^l),fJ.l,fJ.l) ^ f dz z""-' W!,{b,M,s = M'/ z; as{f^l) , 1^% f^l) • (11) 
Jo 

^The subscript b, which labels the parton flavour, should not be confused with the impact parameter b. 
^Throughout the paper, the A^-moments hj^ of any function h{z) of the variable z are defined as = 
dz z^-i h{z) . 
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The resummation structure of ^ can indeed be organized in exponential form, as discussed 
below. 

In the following of this subsection, the subscripts denoting the flavour indices are understood. 
More precisely, we present the resummation formulae in a simplified form, which is valid when 
there is a single species of partons. This simplified form illustrates more clearly the key structure 
and the main features of the resummed partonic cross section. The generalization to considering 
more species of partons does not require any further conceptual steps: it just involves algebraic 
complications, which are discussed in Sect. 12.31 and in Appendix 1X1 

The logarithmic terms embodied in ^ are due to final-state radiation of partons that are 
soft and/or coUinear to the incoming partons. Their all-order resummation can be organized [32] 
in close analogy to the cases of soft-gluon resummed calculations for hadronic event shapes in 
hard-scattering processes [SlHSZl and for threshold contributions to hadronic cross sections [SHI 
IHH] . We write 

WUb, M- as(/i?j), /iL l^l) = K {M, asifxl); M'/^il, MV/i|, M'/Q') 

X expig^iasifil), L- M^fil, M'/Q')} . (12) 

The function does not depend on the impact parameter b and, therefore, it contains all the 
perturbative terms that behave as constants in the limit b —>■ oo. The function Q includes the 
complete dependence on b and, in particular, it contains all the terms that order-by-order in as 
are logarithmically divergent when b oo. This factorization between constant and logarithmic 
terms involves some degree of arbitrariness since the argument of the large logarithms can 
always be rescaled as InM^fe^ = In Q'^b^ + InM^/Q^, provided that Q is independent of b and 
that InM^/Q^ = 0{1) when bM ^ 1. To parametrize this arbitrariness, on the right-hand side 
of Eq. ([T^ we have introduced the scale Q, such that Q ~ M, and we have defined the large 
logarithmic expansion parameter, L, as 

L^ln^, (13) 

where the coefficient b^ = 2e~""^ (7^ = 0.5772 ... is the Euler number) has a kinematical origin 
(the use of 60 in Eq- |T3|) in purely conventional: it simplifies the algebraic expression of Q). 

The role played by the auxiliary scale Q (which we name the 'resummation scale') in the context 
of the resummation program is analogous to the role played by the renormalization (factorization) 
scale in the context of renormalization (factorization). Although the resummed cross section 
does not depend on Q when evaluated to all perturbative orders, its explicit dependence 
on Q appears when is computed by truncation of the resummed expression at some level of 
logarithmic accuracy (see below). As in the case of hr and fip, we should set Q at the central 
value Q = M; variations of the resummation scale Q around this central value can then be used to 
estimate the uncertainty from yet uncalculated logarithmic corrections at higher orders. Note that 
the resummation scale dependence of should not be confused with the 'resummation scheme' 
dependence considered in Ref. jHS]- In fact, as shown in Sect. 12.31 is exactly independent of 
the resummation scheme. 

All the large logarithmic terms a^L"^ with 1 < m < 2n are included in the form factor 
exp{^}. More importantly, all the logarithmic contributions to Q with n + 2 < m < 2n are 
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vanishing. This property, which is called exponentiation, follows pBHHU] from the perturbative 
dynamics of (abelian and non-abelian) gauge theories and from kinematics factorization in impact 
parameter space. Thus, the exponent Q can systematically be expanded as 

g^{as,L;Myfil,MyQ') = L g^'\asL) + {asL; / fi^, Q^) 

+ ^gf{asL;Myfil,M'/Q') (14) 

+00 „ 
• - ■ n—2 



vr 

n=4 



where as = c>^s(/^/j) and the functions g^'^^a^L) are defined such that g^"^"* = when as^^ = 0. 
Thus the term Lg^^^ collects the LL contributions as-Z^'^"'"^; the function g^'^^ resums the NLL 
contributions as-Z^"; g^^^ controls the NNLL terms as-Z^*^"^, and so forth. Note that in the context 
of the resummation approach, the parameter as-Z^ is formally considered as being of order unity. 
Thus, the ratio of two successive terms in the expansion (fT^ is formally of 0(as) (with no L 
enhancement). In this respect, the resummed logarithmic expansion in Eq. (fT^ is as systematic 
as any customary fixed-order expansion in powers of as. 

The function Tiff in Eq. (fT^ does not contain large logarithmic terms to be resummed. It can 
be expanded in powers of as = as(yu|j) as 

H^{M, as; MV/x^, MV/i|, M^Q^) = aP{as, M) [l + ^ 7il^^\M^ / MV/i|, M^/Q^) 

L vr 

+ K^'\M'/^^%M'/^^l,M'/Q') (15) 



n=3 



where a'P = ag a^p'^^ is the lowest-order partonic cross section for the hard-scattering process in 
Eq. (P). 

Two other general aspects of the resummed partonic cross section Wj^ are the factorization 
scale (and scheme) dependence and the process dependence. As discussed below, the form factor 
exp{^} does not depend on both the factorization scale (and scheme) and the specific hard- 
scattering process. 

The hadronic cross section on the left-hand side of Eq. Q is a physical observable and cannot 
depend on the factorization scale fip- In practice, the evaluation of the right-hand side at a certain 
perturbative accuracy introduces the fip dependence of the partonic cross section dapab- This 
dependence is perturbatively balanced by the dependence of the parton densities fa/h{.x,iJ?p). 
Note that the parton densities in Eq. Q do not depend on the transverse momentum qp (or 
on the impact parameter h). Recall also that we implement transverse- momentum resummation 
at the level of the partonic cross section dap at, by using Eqs. © and (fT^ . Therefore, any fip 
dependence of the parton densities cannot introduce any logarithmic dependence on b in the form 
factor exp{^}. In other words, the perturbative expansion (fT3j) of the function depends on 
/Xi?, while the exponent Q of the form factor and its corresponding logarithmic functions (7^'' in 
Eq. (fT^ do not depend on np and on the factorization scheme used to define the parton densities. 
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As explicitly shown in Sect. 12.31 the form factor exp{^} in Eq. (fT^ does not depend on 
the final-state system F produced in the hard-scattering process of Eq. (^. The form factor is 
process independent: it is produced by universal soft and collinear radiation from the QCD partons 
entering the hard-scattering process (when the simplification of considering a single parton species 
is removed, there are various process- independent form factors for the various partonic channels). 
The dependence on the process is fully taken into account by the hard-scattering function 
which embodies contributions produced by virtual corrections at transverse-momentum scales 
qrp M. 

The truncation [VV^j^^ of the resummed cross section at a given logarithmic accuracy is 
defined as follows. At LL accuracy, we include the function g^^^ in the exponent Q and we 
approximate by the Born cross section ap\ At NLL accuracy, we include the functions g^^^ 
and g^^^ and the coefficient Ti.^^^\ At NNLL accuracy, we also include g^^^ and H^^'^\ The 
reason for including both and g)^ at NLL accuracy is that the combined effect of a^Hj^ 

and Lg^^\asL) leads to logarithmic contributions, Ogi^", that are of the same order as those in 

(2) (3) F (1') 

9n i^sL). An analogous observation applies to the inclusion of both g)^ and at NNLL 

accuracy. 

The logarithmic truncation of the resummed component of the cross section can then be 
combined, as in Eq. (jSJ, with the fixed-order expansion of the finite component in Eq. ©. The 
NLL-f-LO result is obtained by supplementing NLL resummation with the LO expansion" at large 
qt- The NNLL-I-NLO result combines NNLL resummation with the NLO expansion at large 
qT- This procedure for combining the resummed and fixed-order approaches exactly satisfies the 
matching conditions in Eqs. (jlj) and ((Zj). Note that the fulfilment of the matching conditions is 
not completely trivial. For instance, if Ti.^^^^ was not included in da^'^^''^ at NLL accuracy, the 
matching condition in Eq. (jZj) would be violated at LO (in other words, Eq. would be violated 
since the qt integral of ['^<5'Faf,'']Lo "^^^^d lead to a non- vanishing finite value when Qt —* 0). 

To reduce the impact of unjustified resummed logarithms in the large-g^ region, we use a 
procedure inspired by that introduced in Ref. [33] to deal with kinematical constraints when 
performing soft-gluon resummation in e~^e~ event shapes. We consider the exponent Q{as,L) of 
the form factor in Eqs. ()12|1 and ()14j) and we perform the replacement 

gias,L) GiasJ) . (16) 

In other words, in the argument of ^(as, L) we replace the logarithmic variable L with the variable 
L defined as 

L.ln(^|!.l) . (17) 

Comparing the definitions in Eqs. (|13p and (|17|). we see that in the resummation region Qb ^ 1 
we have L = L + 0{l/{QbY), and thus the replacement in Eq. (fT^ is fully legitimate** to 
arbitrary logarithmic accuracy. Although the variables L and L are equivalent to organize the 

II We recall that there is a mismatch of notation between the qr distribution at qr ^ M and the total cross 
section. The LO (NLO) term of the finite component of the qr distribution contributes to the total cross section 
at NLO (NNLO). 

**Note that the replacement in Eq. Hl()|l introduces an explicit dependence of da'^p^^'^ on the resummation scale Q. 
Owing to the matching procedure in Eq. (j^Jl, this dependence is balanced by the Q dependence of the da^'^'-*. 



9 



resummation formalism in the region Qb ^ 1, they lead to a different behaviour of the form factor 
at small values of b (i.e. large values of qt)'- when Qb <C 1, we have L — > and exp{^(as, L)} —>■ 1. 
Therefore, performing the replacement in Eq. ()16|) . we reduce the effect produced by the resummed 
contributions in the small-6 region, where the use of the large-6 resummation approach is not 
justified. 

In particular, since exp{^(Q;s5 L)} = 1 at b = 0, using Eqs. © and (fT^ we obtain the relation 



da 



(res.) 



M 



dq, 



which simply follows from the fact that the value at 6 = of the (6-space) Fourier transformation 
of the qT distribution is equal to the integral over qx of the qr distribution itself. Since the 
hard cross section Ti^ is evaluated in fixed-order perturbation theory, the relation ()18|) implies 
that the replacement in Eq. also allows us to implement the perturbative constraint (jH)) on 
the total cross section. More precisely, the integral over qx of the qx distribution dap/dqx at 
NLL+LO (NNLL+NLO) accuracy exactly reproduces the calculation of the total cross section at 
NLO (NNLO). 



The purpose of the transverse-momentum resummation program |26HgO] is to explicitly eval- 
uate the logarithmic functions g^j!^^ of Eq. ()14p in terms of few coefficients that are perturbatively 
computable. As illustrated in Sect. 12.31 this goal is achieved by showing that the all-order resum- 
mation formula (jl4j) has the following integral representation: 



gNiasif^l),L;Myf^l,MyQ' 



' df 



A{as{q^)) In — + S^(as(g')) 



(19) 



where A{as) and -Bat (as) are perturbative functions 



n=4 



(20) 

(21) 



n=3 



The coefficients A^"^ and S^'' are related to the customary coefficients of the Sudakov form factors 
and of the parton anomalous dimensions. This relation is discussed in Sect. 12.31 

Using Eq. ©, the resummed component da^"^^'^ /dq^ of the distribution is fully determined 
by the functions Hf^ and in Eq. (fT2|) . These functions are in turn specified by the perturbative 
coefficients hJ^"^ (see Eq. (US))), and B^j!^^ (see Eqs. ()19|) - (PT|) ). which can be extracted from 
the logarithmic terms in the perturbative expansion of the qr distribution at the n-th relative 
order in as. Therefore, the customary fixed-order computation of the distribution is sufficient 
to obtain the full information that is necessary to explicitly perform resummation at the required 
logarithmic accuracy. 

By inspection of the integration in Eq. (fT^. it is evident that the exponent of the process- 
independent form factor in Eq. (jl2p has the logarithmic structure of Eq. (jl4p . The functions (7^^ 
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depend on the coefficients in Eqs. and (j^H), and the functional dependence is completely 
specified by Eq. (fTT?|l . More precisely (see Eqs. (f^ - P^ ). the LL function g^j^^ depends on A^^\ 
the NLL function g^^^ depends also on A^^-* and B^^\ the NNLL function g^^^ depends also on 
A^^^ and B^^\ and so forth. Starting from the integral representation in Eq. ()19|). the explicit 
functional form of the functions g^j^^ (for arbitrary values of n) can easily be computed by using 
the method that is described in Appendix C of Ref. [T7j . 



The LL, NLL and NNLL functions g^^\ g^^^ and g^^' have the following explicit expressions^^: 



(3) 



g^'\asL) 



[2)1 J M2 M2 



A« A + ln(l- A) 



B 



l3o A 
(1) 



(22) 



N 



+ 



+ 



l3o 
AW 



ln(l - A) 
A 



1-A 



M2 M2 



/5o 

A(3) 



1-A 



-ln2(l 
2 ^ 

A2 



In(l-A)') ln% 
ln(l - A) 



ln(l - A) 



A) + 

■(2) 



A 



B 



N 



1 

A 



+ 



Wl (1 - A)2 
A 



+ In ^ ( 5 



1-A 
(1) A 

'N 



+ 



/5o 1-A 

ln(l - A) 
1-A 



+ 



A 1-A, 

/A(3A 

AW 



(23) 



-A 
2A 



+ 



A^ 



/?o (1-A)^ 



2 (1-A)^ 



2(1- A)2 



;i-2A) In(l-A) 



A 



ln2(l-A) + ln(l-A) 



A 



2(3^(1 - A)2 



- A) 



1-A 

P0P2 - 131 
PI 



+ 



2A 



+ 



ln(l - A) 



0^(1 



A) 



where 



B 



A 



in) 
N 



TT 



5 



/\/f2 

+ /I W In _ 



(n) 



and Pn are the coefficients of the QCD /3 function: 

din as(/i^'' 



(iln/i^ 



/3(as(fi")) 



n=0 



n+1 



(24) 

(25) 
(26) 

(27) 



''^'''Note that the functional form of the functions gjy •* is exactly the same as that of the functions that appear in 
the calculation of the energy-energy correlation in e+e^ annihilation |60j . 
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The explicit expression of the first three coefficients, Pq, Pi and P2, is jEIl 

where Nf is the number of QCD massless flavours and the SU{Nc) colour factors are Ca = Nc 
and Cf = {Nl - 1)/{2N^). 

Note that the functions g^^\asL) in Eqs. are singular at the point A = 1, which in 

terms of the impact parameter corresponds to the value h'^ = h\ = (6q/Q^) exp{7r/(/3oQ!s(/^ij))} 
(i.e. 6l ~ 1/Aqcd, where Aqcd is the momentum scale of the Landau pole in QCD). These 
singularities, which are related (see Eq. (fT^ when b ~ 1/Aqcd) to the divergent behaviour of the 
perturbative running coupling as{q'^)/TT ~ [/5o ln(g^/AQ^^)]~^ near the Landau pole, signal the 
onset of non-perturbative phenomena at very large values of b or, equivalently, in the region of 
very small transverse momenta. 

This type of singularities* is a common feature of all-order resummation formulae of soft-gluon 
contributions. Within a perturbative framework, these singularities have to be regularized. A 
possible regularization procedure consists in introducing a 'minimal prescription', such as those 
introduced in Refs. [32] (in the case of threshold resummation) and [^21 E] (in the case of 6-space 
or joint resummation). In the case of 6-space resummation, other procedures are to use the '6* 
prescription' of Ref. by freezing the integration over b below a fixed upper limit, or more 
simply, to introduce a cut-off at a very large (but smaller than bi) value of b Admittedly, 
when the non-perturbative contributions are sizeable, they have to be properly included, according 
to the prescription used to regularize the singularities. 



2.3 Sudakov form factor, universal form factor and perturbative co- 
efficients 

The 6-space resummation approach was fully formalized by Collins, Soper and Sterman ^Hl 122] 
in terms of perturbative coefficients. Considering the generic hard-scattering process in Eq. (^, 
the transverse-momentum differential cross section in Eq. Q is written as 

^{qr, M,s) = — r db I Mbqr) W^'ib, M, .) + ... , (29) 

where the dots on the right-hand side stand for terms that are not logarithmically enhanced at 
small qT (large b\ Note that Eq. regards the hadronic cross section (and not the partonic 
cross section in Eq. ()in|l ). Therefore, the 6-space function {b^M^s)^ which embodies the 
logarithmically-enhanced terms, depends on the parton densities of the colliding hadrons. The 
all-order resummation of the large logarithms ln(M^6^) in the region Mb ^ 1 is accomplished by 

*Note that these singularities are not related to the presence of factorially-growing coefhcients, such as those 
due to renormalon singularities, at very high perturbative orders. A concise discussion on this point can be found 
in Sect. 3.1 of Rcf. [SHI, in the related context of threshold resummation. 
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showing that the iV-moments W^lb, M) of W{b, M, s) with respect to z = M^/s at fixed M can 
be recast in the following form 



X Y.^-<^^^(''Mb'))C,t,N{as{bl/b')) fa/H,,N{bl/b') h/^^^ , (30) 

a,b 

where fa/h,N{p-'^) are the A^-moments of the parton density fa/h{z, fi"^), and 0"^°^^;. is the lowest- 
order cross section for the partonic subprocess c + c — > F. The function Sc{M, b) is the Sudakov 
form factor of the quark (c = q, q) or of the gluon (c = (7), and it has the following expression"!" : 



5,(M,6)=exp j- r % 



The functions A, i?, C and in Eqs. ()30|) and (jHT|) are perturbative series in as: 



/W"2 

A,(as(g')) In— + i?,(as(g')) 



(31) 



00 

n=l 

00 

^'=("s) = ' (^^^ 

n=l 

00 

a,(as,^) = 5..5(l-^) + ^(^)"cl:)(^) , (34) 



n=l 
n 



i^f(«s) = i + • ^^^^ 

?i=i 

The functions A^, Be and Cah are process independent, while depends on the specific hard- 
scattering process. 

The resummation formulae ()3U|) and ()31|) are invariant under the following 'resummation 
scheme' transformations [S^ : 

i/f(as) - i/f(as) [/^(«s)]"' , 

E>s) - S.(as) - /3(«s) ^^)^^ , (36) 

a mas 

Ca6(as,2) ^ Cab(as,^) [/i(as)]^''^ • 

The invariance can easily be proven by using the following renormalization-group identity (see 
Eq. (EH)): 



|- r'^/3(«s(g^)) 



KaMm = HasiM^)) exp <! - / ^ /3(as(g^)) ^ , (37) 



which is valid for any perturbative function h{as)- 

tin Rcf. 122 the upper limit of the integral in Eq. I|31|l is set to C2M^, where C2 is an arbitrary factor. The 
scale C2Af^ is thus related to the resummation scale in Eq. p9|l . 
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The physical origin of the resummation scheme invariance of Eq. (j3(J|) is discussed in Ref. 
The invariance imphes that the factors , Sc (more precisely, the function Be) and Cab are 
not unambiguously computable order by order in perturbation theory. In other words, these 
factors can be unambiguously defined only by choosing a 'resummation scheme'. The choice of a 
resummation scheme amounts to defining (or Cab) for a single process. More precisely, has 
to be defined for two processes: one process that is controlled, at the lowest perturbative order, 
by qq annihilation (c = q, q) and another process that is controlled by gg fusion {c = g). Having 
done that, the pro cess- dependent factor and the universal (process-independent) factors Sc 
and Cab are unambiguously determined for any other process of the type in Eq. ([T)). 

Note that Eq. (j3Up is usually presented in a form where {as) = 1. Such a form is certainly 
consistent since, by choosing h{as) = (ag) and using the invariance under the transformation 
in Eq. (jHUj) . it is always possible to set H[{as) = 1 on a pro cess- dependent basis. Note that this 
procedure does not correspond to the definition of a resummation scheme. Indeed, the correspond- 
ing Sudakov form factor and the functions turn out to be pro cess- dependent quantities, as 
pointed out by the explicit and general calculation of -Bp'* and C^J{z) in Ref. jnEj- For example, 
in the case of gg fusion processes, the Sudakov form factors for the production of a scalar and a 
pseudoscalar Higgs boson turn out to be different and to have even a different dependence on the 
mass of the top quark. 

Comparing the partonic and the hadronic cross sections in Eqs. (fTUI) and (IHUj) . we see that the 
resummed factors and W^{b, M) are related by 

W^{b,M) = Y,y^Ib,NiP,M■as{^Jil),^ll,^ll) fa/hMP^l) fb/h,,N{f^l) • (38) 

a,b 

To express the resummed partonic cross section in terms of the perturbative coefficients in 
Eqs. we have to use Eq. (jHUI) and substitute the parton densities fa/h,N{bl/b'^) for the 

same parton densities evaluated at the factorization scale fip. The substitution can be done by 
using 

fa/h,N{fJ''^) = ^ Uab,N{fJ'^,fJ'l) fb/h,N{fJ'l) , (39) 
b 

where the QCD evolution operator f/ab^Ar(/i^, /Xq) fulfils the evolution equations 

d Uab,N{f^'^, fJ'o) 



(iln/i^ 



X]^«c,7v(as(/^^)) ^^cfe,iv(/i^/^o) ' (40) 



and 7a6,Af(as) are the parton anomalous dimensions or, more precisely, the A^-moments of the 
customary Altarelli-Parisi splitting functions Pab{o:s, z) 



lab^Nias) = I dz z''-' Pabias,z) = V ( -)"7S;v • (41) 



„l CO 

/ dzz'^-'Pabias,z) = J2{- 

n=l 

We finally obtain [Hn] 

Wi;;,^(6,M;as(/i?j),/xL/^l) = $^a2?p(«s(M2),M)i7f(as(M2))5e(M,6) 



X 5^C,„,,^(as(&o/&'))^ai,^(«s(&o/&')) (42) 

ai, bi 

X Ua^a,Nipl/l?,li\)Ui,^b,Nipl/h^,lil) , 
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which relates the resummed partonic cross section in Eq. (|1U|) to the perturbative coefficients in 
Eqs. fl32|) - (jn^ and the anomalous dimensions coefficients in Eq. ()41|) . 



In the following we explicitly show how Eq. ()42|) is related to the exponential structure of 
Eq. ()12|) in the case with a single species of partons. The general case with partons of different 
flavours is discussed in Appendix El Here we only anticipate that the generalization of Eq. (fT^ 
to the multiflavour case-'- simply involves a sum of exponential terms, namely 

K,,Ab,M;as{i^l),i^%f^l) = 5^^;;?/ (M,as(/x?i);MV/iL^V/i|,MVg') 

{1} 

xexp{g{j},r^{as{fil),L;My^l,MyQ^)} , (43) 
where the index {/} labels a set of flavour indices (which is precisely specifled in Appendix 

Within the simplifled treatment in which there is a single species of partons, the resummed 
partonic cross section in Eq. (j^^ can easily be recast in the factorized exponential form of Eqs. (fT^ 
and (fnH) . To this aim, we flrst use the identity (pTTjl with h{as) = CN{as) to replace CAr(as(&o/^^)) 
in Eq. (jl^ in terms of Cat (as (M^)). Then, we insert in Eq. ()42|) the solution of the evolution 
equation (jlUJ: 

UM/b^l^l) = exp I - r % 77v(«s(g')) 1 . (44) 

We flnally obtain the exponential form in Eq. ()19j) . where the perturbative function A{a^) is 
exactly the perturbative function in Eq. ()32j) . and the function i?Ar(as) is given as follows in terms 
of the perturbative functions in Eqs. (jTTj) . (jH^ . (jMjl and (jUJ: 

Sat as = 5 as + 2/5 as + 27^(as) . (45) 

a mas 

The expression of the hard-process function in Eq. (fT^ is 

H^(M,as(/x^);MV/iL^V/i|,MVg') = aP{a^{M'),M) H^{as{M^)) C],{as{M')) 



X exp 



dq 



2 



A(as(g')) In — + S^(as(g')) 



Note that, as discussed in Sect. 12. 2t the form factor exp{^} and, hence, the perturbative 
functions v4(as) and Bjq{a^) in Eq. (fT^ do not depend on the factorization scale ^ip- As a 
consequence, the functions A{as) and -Bat (as) are also independent of the factorization scheme 
used to deflne the parton densities. Since, as is well known, the anomalous dimensions 7af),7v(as) do 
depend on the factorization scheme, the relation implies that both the perturbative functions 
Bc[aQ) and Cab(as) depend on the factorization scheme in such a way that i?Ar(as) turns out to 
be factorization-scheme independent. 

As anticipated in Sect. 12.21 the form factor exp{^} does not depend on the flnal-state system F 
produced in the hard-scattering process. From Eqs. (|Tn|l and pHjl . this independence is a simple 

■fin the multiflavour case, Eq. H12f) directly applies to the flavour non-singlet components of the resummed 
partonic cross section. 
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consequence of the process independence of each of the perturbative functions Ac{as), Bc{as), 
Cabins) and 7a6,Af(as)- 

The relation ()45|) also implies that the form factor exp{^} does not depend on the resumma- 
tion scheme used to express the various factors in the resummation formulae ()3U|) and ()31|) (we 
recall that the customary Sudakov form factor Sc{M, h) in Eq. (p?T|l does instead depend on the 
resummation scheme). It is indeed straightforward to show that the function i?Ar(as) in Eq. 
is invariant under the resummation-scheme transformations in Eq. ()36|) . 

Unlike the form factor exp{^}, the non- logarithmic function Ti^ in Eq. pUj) explicitly depends 
on the factorization scale ^p, on the factorization scheme (through Caf,,Ar(as) and 7afe,Ar(as)) and 
on the final-state system F (through af) and H^). Nonetheless, Ti^ does not depend on the 
resummation scheme, since the factor {a^)Clf{a^) is invariant under the transformations in 
Eq. (P^. 

The universal (i.e. independent of the process and of the factorization and resummation 
schemes) perturbative function Ac (as) in Eqs. PU)) and (jH^ is known up to 0{a1). The LL and 
NLL coefficients A'i^ and Af^ are |Si El 

where Cc = Ci? if c = g, g and Cc = Ca if c = g. The NNLL coefficient Af is not yet known. In 
our quantitative study of transverse-momentum resummation at NNLL accuracy (see Sect. El), we 
assume that the value of Ac is the same as the one [HSlinS] that appears in resummed calculations 
of soft-gluon contributions near partonic threshold. This assumption is based on the fact that the 
two coefficients A^c^ and A^c^ in Eq. (j47|) are exactly equal to those of the related perturbative 
function that controls threshold resummation j^H] in the MS factorization scheme. Note, however, 
that the two soft-gluon functions Ac{as) do not necessarily coincide at high perturbative orders 
since, for instance, the soft-gluon function for transverse-momentum resummation is universal 
while the soft-gluon function for threshold resummation depends on the factorization scheme. 

The first-order coefficient -B^^]^ of the universal perturbative function B]\f{as) in Eqs. (j2II) and 
(HSl) is 

= Bi'^ + 27£:^ , (48) 

where ^ 03 

= = -lCj. , i?« = (IIC^ - 2Nf) . (49) 

Note that, since the LO anomalous dimensions 7£'';v are universal, the NLL coefficients B^^^ in 
Eq. (j49p are themselves independent of the factorization and resummation schemes. 

The universal second-order coefficient B'^]^ in Eq. (^3]) is 

55 = 5p)-2/3oC«, + 27S, , (50) 
or, equivalently, by performing the inverse Mellin transformation to z-space: 

Bi'\z) = 6{1 - z) - 2(3, C^^{z) + 2P|)(z) . (51) 



67 vrA ^ 5^^ 

\Ca - -Nf 

18 6 9 ^ 



(47) 
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The value of the quark coefficient Bg can be obtained by using the results of Ref. jHZI for 
the coefficients Bg^^ and Cgl\z) of the DY process. These results are conffimed by the general 
(process-independent) calculation of Ref. [201 , which considers both the gg-annihilation and the 
gluon fusion channels. From the results of Ref. [HE] we obtain the value of the gluon coefficient 
B'g^\ and we can also explicitly check the universality of both Bg and Bg . To write down the 

'^(2) (2) 

expression of i?c , we recall that the second-order term Pec (z) of the Altarelli-Parisi splitting 
functions Pcc{(Xs,z) has the following general dependence on z: 

P^{z) = + 6{l - z) \ 7f + i^f , (52) 

where Ac is the coefficient in Eq. (^Tj), 1/(1 — z)+ is the customary 'plus'-distribution and 
pj2)regJ^^^J denotes all the remaining and less singular (when z ^ 1) contributions to Pcc\z). The 
explicit expressions of Pc?'^'^^{z) and of the constants "yi^^ can be found in the literature (see e.g. 
Ref. ini). Using the notation of Eq. the universal NNLL coefficient bP is ^ 

B?\z) = + 5(1 - z) P.Cc'^ + 2P(2)reg(^) + 2/3oAU^) , (53) 

{l-z)+ 6 

where 

nv^) = -\ (1 - ^) ' ^.y^) = • (54) 

The ffist-order coefficients C'^g and C^q in Eq. (jH^ do not depend on the process and on 
the resummation scheme, and were ffist computed in Refs. [HZl and [33], respectively. Their 
expressions in the MS factorization scheme are 

^SH^) = Cf^{z) = \z{l-z) , C^/J{z) = &S{z) = \Crz . (55) 



The flavour-diagonal ffist-order coefficients C'^q and C^g and the coefficients Hq and Hg '•^^ 
depend on the resummation scheme. The dependence on the resummation scheme is cancelled 
in the perturbative coefficients of the hard-process function Tifj. For example, by expanding 
Eq. (j46j) in powers of as(yUif,)) we obtain the following expression for the ffist-order coefficient 
KJ^^^ of Eq. (dSD: 




i? = In — , iF = ln— , £q = In — . (57) 



where we have defined 

— , = In — , £q = In — 
/ijj Hp LJ 

The coefficient 'H^^^^ depends on the process and is explicitly known for several processes (see 
Ref. [HH] and references therein). 

To complete the resummation program at NNLL, the coefficient 7i^*^^^ is also needed. This 
coefficient is not known in analytic form for any hard-scattering process. Nonetheless, within our 
resummation formalism, it can be determined for any hard-scattering process whose corresponding 
total cross section is known at NNLO. This point is discussed in detail at the end of Sect. 12.41 
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2.4 The finite component 



The finite component da^p'^jj / dq'^ of tlie transverse- momentum cross section is computed at a given 
fixed order in as according to Eq. ©. To implement Eq. (0), we liave to subtract ['^<^Faf)'']f o 
from [d^Fabli^ - 

As discussed in Sects. 12.11 and 12. 2| tlie finite component daf'^^/dq^ does not contain any 
perturbative contributions proportional to 6{q^) (these contributions and all the logarithmically- 
enhanced terms at small qx are included in da^p''^^ / dq^) . Therefore, when computing ['^'3"i?ab'*]£Q 
according to the subtraction procedure in Eq. (jHl), we can consistently neglect any terms propor- 
tional to S{q^) both in [(ia^abjf ^ and in ['^o'^afe^lf o ■ "^^^^ formally equivalent to the evaluation 
of both [(iciFafejfQ and ['^^"^afe'*] f o large-gr region (or, more precisely, in the region where 

gr 7^ 0). The expansions of ['^0"Fa6^]f o ^^^^ ^^"^ second perturbative order thus give 



da 



(fin.) 
Fab 



da 



T 
(fin.). 



- dq, 



LO 



da 



Fab 



dq. 



J NLO 



dq^ 
da 



J LO 



da 



(res.) 
Fab 



Fab 



dqp J NLO 



dq; 
da 



T 

(res.) 
Fab 



LO 



dq. 



J NLO 



(58) 
(59) 



where the subscript LO (NLO) denotes the perturbative truncation of the various cross sections 
at the leading order (next-to-leading order) in the region where qx 0. The extension of Eqs. (j^Hjl 
and ()59|) at still higher perturbative order is straightforward. 

The contributions [c/ai^afejf^ on the right-hand side of Eqs. (fSHj) and are obtained by 
computing the customary perturbative series for the partonic cross section at a given fixed order 
(f.o.=LO, NLO, ...) in ag. The fixed-order truncation ['^o'^ab''] f ^ of the resummed component 

is obtained by perturbatively expanding the resummed component da^^^^'^ in Eqs. (fTIH) . To this 
purpose, we define the perturbative coefficients S*^"^ as follows: 



Wf,(6, M, s- as, fil, fil, Q') = ^S!f(«s, M) { 6,a 5,, 5(1 

c 

oo 



n=l 



,F(n) 



m2 m2 



cc^ab \ ^' ,,2 ' ,,2 ' 



^-R fJ'p 



where z = M^/s, 



as 



«s(/i?j) 



a 



(0) 
cc, F 



'as,M) 



a: 



^cc F (^) ^'^'i' ™- general, the power p^F 



depends on the lowest-order partonic subprocess c + c — > F. In Eq. (I6(J|) . is the resummed 
cross section on the right-hand side of Eq. (jTU)). Note, however, that Eq. (jUUI) depends on the 
resummation scale Q^. The dependence on the resummation scale has been introduced in Eqs. fjlOj) 
and (fT^ through the replacement in Eqs. (fTBj) . The perturbative coefficient S*^") on the right-hand 
side of Eq. (jUn|) is a polynomial of degree 2n in the logarithmic variable L defined in Eq. ()17|) . 
The coefficients S'-"^ vanish by definition when L = (i.e. when 6 = 0), and the 6-independent 
part of ^(6, M) is embodied in the coefficients 'H^'^\ 
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The perturbative expansion of Eq. (fT^ or, more precisely, of Eq. (li^ gives 



''cc<—ab 



''cc^ab 



(61) 



(62) 



where the dependence on the scale ratios M^/ fij^, M^/ fip and M^/ is understood. The extension 
of Eqs. and to the higher order terms E^^^f,(2;,L) with n > 3, is straightforward. The 
6-independent coefficients ^^'''\z),n^ ^^\z),E^ ^^'''\z) and n^^^\z) are more easily presented 
by considering their A^-moments with respect to the variable z. We have 



.F(l;2) 
-'cc<— afe, N 



(63) 



''cc<— ab, N 



1 ca, N 



(64) 



n 



F(i) I M2 



+ + S-cbC^],, + [S^al^N + S-cdl]N) i^F ' ^q) , (65) 



(66) 



-'cc*—ab, N 



{M'/Q') 



(67) 



p(2;2) M2 



cc^ab,N \ 2 ' ,,2 ' /02 



/^F 



- - ^i'^ 

2 

1 



F(l) 



M2 M2 M2 



cc^afe, I "X' T;2"' 7^ ) ~ ScaScb {^R " ^q) 



(68) 



^(2;l) fM^ M2 M2\ _ ^(1.1) 



''cc'—ab 



ai,6i 



M2 M2 M2 



(69) 



(2) 

ca, Af 
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(1) 

cb,N 



F{2) (M^ _ , r o-F(2) , c ^(2) c ^(2) ^(1) ^ 

'^cc^ab,N\ 2 ' 2 ' n2 i — '■> ca<J cb c '■'ca'-^ cb. N ^ ca, N ^ ca, N t 



cc'^ab,N\ 2 ' ,,2 



M2 M2 M2 

^5 0~? ^9 P'O*- 



"T r, 2-^ ''-cc^aibi, Af I , 2 ', 2 ' n2 i "I" "^caj Ocbi "I" "^cai "-^cbi , AT "cfei '-^cai 



- Search PcF [ ^Po^R + Pi 



R 



(70) 



The right-hand side of Eqs. ()63 p - (f7m) is expressed in terms of the resummation-scheme independent 
coefficients given in Sect. 12.31 and of the logarithms iqjip and in defined in Eq. (jSTj). To exphcitly 
exhibit the independence of the resummation scheme we can, for example, rewrite the contribution 
in the third line of Eq. (|69|) in terms of the resummation-scheme independent coefficients B'^^ (see 
Eq. and C^J^^ with a^b (see Eq. 



+ 5ca (1 - S,,) (7g,V - PoC^J,) + S,b (1 - 5ca) (7g7V ' PoC^l] N 



(2) 

ca,N 



(71) 



Inserting Eqs. (j60|) - (j62j) in Eq. (jlOj) . performing the integral over the impact parameter b, 
and removing the contributions proportional to S{q^) (for example, all the contributions coming 
from H^Jl^li^ in Eq. (plj) ). we obtain the following expressions for the fixed-order contributions 
['^^Fab^]f o right-hand side of Eqs. (j3Hj) and (j^Uj) : 



da 



(res.) 



dq, 



M,s 



M2 



LO 



«s(/i|) ^ 

vr Q2 



E<^^("s(^?j),M) 



X 



/2(g./Q) + sf,a^ ( ^ ) IMQ) 



(72) 



2 J NLO 



dq^ 



da 



(res.) 
Fab 



dq^ 



{qr, M, s; asifin), fin, fip^Q 



LO 



+ 



«s(/^|) 



TT 



On the right-hand side of Eqs. (f?^ and (fTHj) . the dependence on qx is fully embodied in the 
functions Iniqr/Q), which are obtained by the following Bessel transformation: 



In{qT/Q) = Q' / db-Mbqr) In" 



/)2 



(74) 
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The term ln"(l + Q'^b'^/bl) = in the integrand comes from the replacement L ^ L (see 
Eq. (fT^ ). In customary implementations of 6-space resummation, one has to consider the Bessel 
transformation of powers of In^lQ^b"^ /bl) = L", which can be expressed in terms of powers of 
ln"'((5^/g|'). The functions Inilr/Q) have instead a more involved functional dependence on g^. 
As shown in Appendix IHl this functional dependence can be expressed in terms of K^{qT/Q), 
the modified Bessel function of imaginary argument that is defined by the following integral 
representation: 

POO 

KUx)= I rft e-^''=°"^*coshz/t . (75) 



We conclude this section with some observations on the hard-scattering function H^^^^i^. This 



function is resummation-scheme independent, but it depends on the specific hard-scattering sub- 
process c + c ^ F. The coefficients Ti^J^l^i^ of its perturbative expansion can be determined 
by performing a customary perturbative calculation of the qt distribution in the limit qt — > 0. 
Moreover, as discussed in Sect. 12.21 within our resummation formalism Ti.^ controls the strict 
perturbative normalization of the corresponding total cross section (i.e. the integral of the qt 
distribution). This property can be exploited to determine the coefficients Ti-^J^l^fj in a different 
manner, that is, from the perturbative calculation of the total cross section. 



To illustrate this point we consider the total cross section, o"^*^, at the partonic level 



a^°*,(M,s;«s(/iJj),/ii/iF) 



dap 



^^2 (9t,M, s;as(/i^),/i^,/i^) 



(76) 



and we evaluate the qt spectrum on right-hand side according to the decomposition in terms of 
'resummed' and 'finite' components (see Eq. (jS))). Then we use Eq. (fTHj) to integrate the resummed 
component over qt, and we obtain 



"tot 



M2 
s 



T^ab + 



"^^^ dql 



(77) 



This expression is valid order by order in QCD perturbation theory. Once the perturbative 
coefficients of the fixed-order expansions of o"^*;,, TifJ, and da^p^jj /dq^ are all known, the relation 
(f77|) has to be regarded as an identity, which can explicitly be checked. Note, however, that since 
the fixed-order truncation [da^p^b /^^rlf o does not contain any contributions proportional to 
5(g|.), [do-^p^^ /dq^ does not explicitly depend on the coefficient Ti.^J'^^ (see Eqs. (I58|) and (f7^ ^ 



LO 



Analogously, [d&p °^ / dq- 



NLO 



does not explicitly depend on the coefficient H^f^ 



F{2) 



[see Eqs. §^ 



F{n) 
ab 



and (f73j)). and so forth. Therefore, Eq. (f77|) can be used to determine the N'^LO coefficient H 
from the knowledge of o"^*;, at N°LO and of da^p^^ /dq^ at N°~^LO, without the need of explicitly 
computing the small-g^ behaviour of the spectrum dapab/dq^ at N°LO. For example, at NLO 
Eq. (f77j) gives 



TT S ^ 



a 



M^-u^i) (M' M2 M2 M2 



(0) 



a';l,{M,s;as,fiji,fi'p) 



NLO 



LO 



dqp 



da 



(fin.) 
Fab 



dq'l 



{qT,M,s;as,f^n,f^l,Q^ 



(7J 



J LO 
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where as = as(/i|;) and we have used 



a^?l,(M, s; as) =5(1- M^/s) ^ <V(as, M) 4a 5a 

J LO 

c 

At NNLO Eq. dZZj) gives 

/as\2 (0) ,.^^F(2) ( M' M' M'\ 



cr 



Fab 



NNLO 



a 



„tot 
Fab 



NLO 



da 



(fin.) 



NLO 



da 



(fin.) 
Fab 



LO 



and the generahzation at still higher orders n {n > 2) is 



TT 



F(n) 
cc^ab 



-tot 
"Fab 



N°LO 



da 



a 



-tot 
Fafe 



Nn-iLO 



(fin.) 
Fab 



dq^ 



Nn-iLO 



da 



(fin.) ^ 
Fab 



(79) 



Nn-2LO 



^0) 



In our study of the transverse-momentum spectrum of the Higgs boson at NNLL accuracy (see 
Sect. El), we use Eq. (jHUI) to obtain a numerical value for the corresponding perturbative coeffi- 
cient n^'^\ 



3 The qr spectrum of the Higgs boson at the LHC 

In this section we apply the resummation formalism described in Sect. El to the production of the 
SM Higgs boson at the LHC. 

We consider the gluon fusion production mechanism gg H, whose Born level cross section 
in Eqs. 1^ and ^ is 

aSJ^(as, Mh) = 6,, 4 (M^; M,) , (82) 

where Mt and Mf, denote the masses of the top and bottom quark, which circulate in the heavy- 
quark loop that couples to the Higgs boson. In our numerical study we use Mt = 175 GeV 
and Mh = 4.75 GeV. The expression of a^^^Mff; Mt, Mf,) can be found, for instance, in Eq. (3) 
of Ref. [T7j. Though the Born cross section is evaluated exactly, i.e. including its dependence 
on the top- and bottom-quark masses, the computation of the higher-order QCD corrections is 
performed in the framework of the large-M^ approximation. More precisely, we proceed as in 
Ref. ^7]: we first compute dan/dqT in the large- limit and then we rescale the result by the 
factor a^°\MH;Mt,Mb)/a!^\ where is obtained from a^^^Mn; Mt, Mb) by setting Mb = 
and Mt/Mu ^ cxd. As recalled in Sect. ^ this implementation of the large- approximation is 
expected to produce an uncertainty that is smaller than the uncertainties from yet uncalculated 
perturbative terms from higher orders. 
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We compute the Higgs boson differential cross section da/dq^ at the LHC {pp colhsions at 
■s/s = 14 TeV) and present quantitative resuhs at NLL+LO and NNLL+NLO accuracy. 

As discussed in Sect. 12.21 at NLL+LO accuracy the resummed component in Eq. (fT^ is 
evaluated by including the functions (7*^^^ and g^^^ in Eq. (fT^ and the coefficient H^^^^ in Eq. (fTS|) . 
and then it is matched with the fixed-order contribution evaluated at the LO (i.e. at 0{ag)) in 
the large-gr region. The functions g^^^ and g^^^ are process independent and given in terms of the 
universal coefficients A^^\A'''^'> and i?^'* (see Sect. I2.3|l . The flavour off-diagonal part of 'H^gllab.N 
is also process independent and given by Eq. (jHSJ; setting fip = Q = Mh, we simply have 

njH{l) _njH{l) _ ^(1) _ 1 (c,n\ 

'^gg^gq,N - l^gg^qg,N " ^ gq,N " 2(A^ + 1) ^ ' ^ ' 

where the coefficient C^gqi^ is the Mellin transformation of Eq. Of course, these process- 

independent coefficients are exact, i.e. they are not affected by the large-Mj approximation. The 
flavour diagonal coefficient Hggjgg is instead process dependent; therefore it depends on Mt and, 
in the large- approximation, it is given by 



Kg^gg,N = Hf^'' + KJ^n = ^ [(5 + - SC.] = ^ (11 + 3.^) , (84) 

where, for simplicity, the scale-dependent terms have been dropped (i.e. we have set fiR = fip = 
Q = MHm Eq. 

At NNLL-t-NLO accuracy the function g)^ and the coefficient Ti^y have also to be included 
in the resummed component of the qt cross section, and the finite component has to include the 
fixed-order contribution to the cross section evaluated at the NLO (i.e. at 0{ag)) in the large-gr 
region. The process-independent function g^^^ depends on the universal coefficients A^^'^ and i?^^ 
(see Sect. 12. 3|) . The scale-independent part of the coefficient H^^'^^ (its scale-dependent part can 
be obtained from Eq. (j7(Jj) ) is not known in analytic form. We thus exploit Eq. (j8(jp . which follows 

17(0') 

from the constraint of perturbative unitarity, to extract the numerical value of from the 

knowledge of the total cross section at the NNLO [131. The scale-independent part of T^^g'^gg n 
can be written as 

. = H^''' +<1n + (C^'JS + 2<« (85) 
19 1 \ , M2 

where the Mj-dependent contribution on the right-hand side is obtained from the results in 
Refs. [iniinH], and cat does not depend on Mf in the large- approximation. Since from Eq. 



'^gg^gg, n 



we know that C^gg m actually independent of A^, the A^ dependence of cat can only follow from 
that of Cfl ^ in Eq. (j85|). Using Eq. (j80|l and the NNLO total cross section, we find that the 

-F (1'] 

fiavour off-diagonal terms in Ti^^ J^j^ ^ can numerically be neglected, and that the coefficient cat 
in Eq. (jHUj) can numerically be approximated by an A^-independent value, cat ~ 178.75. This 
numerical approximation is pretty good, since the integral of the NNLL+NLO spectrum repro- 
duces the NNLO total cross section to better than 1% accuracy in a wide Higgs mass range, 
100 GeV<Mj^<300 GeV, at the LHC. 
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We recall that the functions gjsi (X) are singular when A — 1 (see Eqs. (I22 |l - (j24|l ). The 
singular behaviour is related to the presence of the Landau pole in the perturbative running of 
the QCD coupling asiq"^)- As mentioned at the end of Sect. \2.2\ a practical implementation of the 
resummation procedure requires a prescription to deal with these singularities. In our numerical 
study we follow Ref. [^2] and deform the integration contour in the complex b space. In particular 
we choose the two integration branches as 

b = {cos (f)±i sin (j))t , t G {0, 00} . (87) 

We have checked that the result is very mildly dependent on the choice of 0. We have also used 
the simpler procedure of integrating over the real 6-axis, using a sharp cut-off at a large value of 
b, and checking the independence of the actual value of the cut-off. We found that the numerical 
differences between the results obtained by these two procedures are negligible. 

Our complete calculation of the qx spectrum of the Higgs boson at the LHC is implemented 
in the numerical code HqT, which can be downloaded from [^2] together with some accompanying 
notes. This code is a slightly modified and numerically improved version of the code used in 
Ref. the most important difference regards the computation of the finite component. In 
Ref. P we used the Monte Carlo program of Ref. [20] to compute the fixed-order contribution 
to the qt cross section at LO and NLO. Here we have implemented the analytic calculation of 
Glosser and Schmidt [22] • Although the two methods are in principle equivalent, the use of the 
analytic calculation allows us to achieve a faster numerical stability in the small-gr region. In the 
next subsection we present a selection of numerical results that can be obtained with our code. 
We also include a discussion of theoretical uncertainties. 

3.1 Numerical results at the LHC 

To compute the hadronic cross section, we use the MRST2004 set [70] of parton distribution 
functions. As for the perturbative order of the parton densities and a^, at variance with Ref. P, 
we adopt here the following choice. At NLL+LO we use NLO parton densities and 2- loop as, 
whereas at NNLL-I-NLO we use NNLO parton densities and 3-loop as- This choice is perfectly 
consistent in the small qt region, since the corresponding partonic cross section is dominated by 
the resummed component evaluated at NLL and NNLL accuracy, respectively. The choice is fully 
justified also at intermediate values of qt, where the calculation of the partonic cross section is 
driven by the small-g^ resummation and strongly constrained by the total cross section at NLO 
and NNLO, respectively. At large values of qt, Qt ~ Mh, our evaluation of the partonic cross 
section is dominated by the fixed-order contributions at LO and NLO, respectively. Therefore, 
our choice introduces a formal mismatch with respect to the customary use of parton densities 
and as- However, as shown and discussed later in this subsection, this formal mismatch does not 
lead to any inconsistencies at the quantitative level. 

The NLL+LO spectrum with = 125 GeV is shown in Fig. In the left-hand side, the 
full NLL+LO result (solid line) is compared with the LO one (dashed line) at the default scales 
I^F = IJ^R = Q = Mh- We see that the LO calculation diverges to +00 as qt 0. The effect 
of the resummation, which is relevant below qr ~ 100 GeV, leads to a physically well-behaved 
distribution: it has a kinematical peak at qj- ~ 12 GeV and vanishes as qx 0. The LO finite 
component of the spectrum (dotted line), which is defined in Eq. (j58|) . is also shown: as expected 
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Figure 1: The qx spectrum at the LHC with Mh = 125 GeV: (left) setting fJ^R = fJ^F = Q = ^h, 
the results at NLL+LO accuracy are compared with the LO spectrum and the finite component of 
the LO spectrum; (right) the uncertainty hand from variations of the scales /i/j and fip at NLL+LO 
accuracy. 

it dominates when qt ~ Mh and vanishes as gr 0. Note, however, that the contribution of the 
finite component is sizeable in the intermediate-g-r region (about 20% at qx ~ 50 GeV) and not yet 
neghgible at small values of qx (about 8% around the peak region). This underlies the importance 
of a careful and consistent matching between the resummed and fixed-order calculations. In 
the right-hand side of Fig. ^ we show the NLL-I-LO band as obtained by varying fip and fiR 
simultaneously and independently in the range 0.5Mh < fJ'F,fJ'R < 2Mh with the constraint 
0.5 < fJ.p/fJ'R < 2 (the resummation scale is kept fixed at Q = Mh)- The scale dependence 
increases from about ±15% at the peak to about ±20% at qx = 100 GeV. The integral over qx 
of the NLL±LO spectrum is in agreement with the value of the NLO total cross section to better 
than 1%, thus proving the numerical accuracy of the code. 

The NNLL±NLO results at the LHC are shown in Fig. |21 In the left-hand side, the full result 
(solid line) is compared with the NLO one (dashed line) at the default scales lip = I^r = Q = Mh- 
The NLO result diverges to —oo as and, at small values of qx-, it has an unphysical peak 

(the top of the peak is above the vertical scale of the plot) that is produced by the numerical 
compensation of negative leading logarithmic and positive subleading logarithmic contributions. 
The resummed result is physically well-behaved at small qx- The NLO finite component of the 
spectrum (dotted line), which is defined in Eq. vanishes smoothly as qx — > 0; its contribution 
amounts to about 10% in the peak region, about 17% at qx ~ 25 GeV and about 35% at qx ~ 
50 GeV. This shows both the quality and the relevance of the matching procedure. 

We find that the contribution of A^^'^ (recall from Sect. I2.3l that we are using an educated guess 
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Figure 2: The qt spectrum at the LHC with Mh = 125 GeV: (left) setting fJ^R = fJ^F = Q = Mh, 
the results at NNLL+NLO accuracy are compared with the NLO spectrum and the finite component 
of the NLO spectrum; (right) the uncertainty band from variations of the scales fipt and fip at 
NNLL+NLO accuracy. 

on the value of the coefficient A^^^) to the resummed component can safely be neglected. The 
coefficient H^^"^^ contributes significantly, and enhances the qt distribution by roughly 20% in the 
region of intermediate and small values of qt- The NNLL resummation effect starts to be visible 
below qx ~ 100 GeV, and it increases the NLO result by about 25% at = 50 GeV. 

The right-hand side of Fig. |21 shows the scale dependence computed as in Fig. [T] The scale 
dependence is now about 8% at the peak and increases to about 20% at qx = 100 GeV. 

To better illustrate the main features of the dependence on the scales fiR and /Xi?, we present 
numerical results at two fixed values of qx in Figs. 01 and El In Fig. El we show our results at 
qT = 50 GeV and Mh = 125 GeV. The scale dependence is analysed by varying the factorization 
and renormalization scales around the default value Mh- The plot on the left corresponds to 
the simultaneous variation of both scales, = f^R = X^h, whereas the plot in the centre (on 
the right) corresponds to the variation of the factorization (renormalization) scale fip = Xf Mh 
{Hr = xr Mh) by fixing the other scale at the default value Mh- 

As expected from the QCD running of as, the cross sections typically decrease when ^r 
increases around the characteristic hard scale Mh, at fixed hf = Mh- In the case of variations 
of fip at fixed fiR = Mh, we observe the opposite behaviour. This is not unexpected, since 
when Mh = 125 GeV the cross section is mainly sensitive to partons with momentum fraction 
X ~ 10~^, and in this x-range scaling violations of the parton densities are (moderately) positive. 
Varying the two scales simultaneously {fiF = (^r) leads to a partial compensation of the two 
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Figure 3: Scale dependence of the LHC cross section for Higgs boson production (Mh = 125 GeV) 
at qt = 50 GeV. Results at a) (upper) LO, NLL+LO and b) (lower) NLO, NNLL+NLO accuracy. 

different behaviours. As a result, the scale dependence is mostly driven by the renormalization 
scale, because the lowest-order contribution to the process is proportional to ag, a (relatively) 
high power of as- 

Comparing the LO with the NLL+LO results and the NLO with the NNLL+NLO results, we 
see that the scale dependence of the resummed results (solid lines) is smaller than that of the 
corresponding fixed-order results (dashed lines): the LO and NLL+LO curves have a comparable 
slope, but the NLL+LO results are higher; the NLO and NNLL+NLO results have smaller differ- 
ences, but the slope of the NNLL+NLO curve is flatter. In summary, resummation reduces the 
scale dependence of the fixed-order calculations also in the region of intermediate values of qt- 

In Fig. m we report analogous results at a smaller value of qt, namely qt = 15 GeV. The 
qualitative behaviour is similar to the one in Fig. |21 In this region of small transverse momenta 
the fixed-order result is no longer reliable (see Figs.[T]and|21), but its relative scale dependence does 
not increase and is even smaller than at qt = 50 GeV. This is due to the fact that the fixed-order 
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Figure 4: Scale dependence of the LHC cross section for Higgs boson production (Mh = 125 GeV) 
at qt = 15 GeV. Results at a) (upper) LO, NLL+LO and b) (lower) NLO, NNLL+NLO accuracy. 

cross section is much larger than at higher values of qt- The slope of the resummed results (solid 
lines) is sizeably flatter than that of the corresponding fixed-order results (dashed lines). We also 
notice a slight reduction in the scale dependence of the resummed results compared to Fig. El 
especially at NNLL+NLO accuracy. 

In Fig. El the NLL+LO and NNLL+NLO bands shown in Figs. Q and |21 are compared. We 
see that the NNLL+NLO band (solid lines) is smaller than the NLL+LO one (dashed lines) and 
overlaps with the latter at qt ^ 100 GeV. This suggests a good convergence of the resummed 
perturbative expansion. This result is confirmed by the inset plot, that shows the NNLL+NLO 
band normalized to the NLL+LO result at central value of the scales. This g^-dependent K factor, 

7^/ N daNNLL+NLoifJ^F, (J'R) 

^ y^T) - -1 7 T-j-T , [oo) 

dcTNLL+LOKl^F = (J'R = Mh) 

is stable, around the values 1.1-1.2, in the central region of the inset plot, and it increases (de- 
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creases) drastically when gj'>50 GeV {qT^2 GeV). In the large-g^ region, the effect of pertur- 
bative higher-order corrections is known to be important |2nH221- At very small values of qt, 
non-perturbative effects are definitely expected to be relevant. We observe that a naive rescaling 
of the NLL+LO result by a constant (i.e. independent of qt) K factor would not reproduce the 
NNLL+NLO result over the entire gr-range. 
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Figure 5: Comparison of the NLL+LO and NNLL+NLO bands (Mh = 125 GeV). The inset plot 
shows the NNLL+NLO band normalized to the central value of the NLL+LO result. 

The nice convergence of the resummed perturbative expansion suggested by Fig. El should be 
contrasted with the results in Fig. IHl where the corresponding fixed-order bands, computed as in 
Fig.El are shown. The results in Fig. IHl have no physical significance in the small-g^- region, owing 
to the non- convergence of the fixed-order expansion herein. When qt ^ 25 GeV, we see that the 
scale dependence of the NLO (LO) result is larger than the one of the corresponding NNLL+NLO 
(NLL+LO) result in Fig. |3| More importantly, we see that the LO and NLO bands do not 
overlap. This implies that the scale dependence enclosed by these bands certainly underestimates 
the true theoretical uncertainty from missing higher-order terms. Equivalently, we can say that 
the uncertainty of these fixed-order calculations is more reliably estimated by performing scale 
variations over a range of scales that is wider than that used in Fig. IHl AH this indicates a poor 
convergence of the fixed-order perturbative expansion at intermediate values of qt 

As mentioned at the beginning of this subsection, in our resummed calculations at NLL+LO 
and NNLL+NLO accuracy we use parton densities and as at perturbative orders that are different 
from those customarily used in fixed-order calculations at LO and NLO, respectively. Indeed, the 
consistent procedure at large values of qr would be to use LO densities with 1-loop as at the 
LO, and NLO densities with 2-loop as at the NLO. We have also explained why our procedure 
is justified in the intermediate-gT region, and we have postponed the discussion on the large-g^ 
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Figure 6: Comparison of the LO and NLO bands (Mh = 125 GeV). 

region. To come back to this point, in Fig. [7|we compare our NLL+LO and NNLL+NLO results 
with the customary NLO results, which are obtained by using NLO parton densities and 2- loop as- 
We also include the corresponding bands, computed from scale variations. In the left-hand side 
we see that in the intermediate-gr region our NLL+LO result catches the bulk of the NLO effect. 
Obviously, at large qt, the inclusion of NLO corrections is necessary. In the right-hand side, the 
calculations at NNLL+NLO accuracy and at the NLO are compared. In spite of the fact that the 
two calculations use different parton densities and as, the corresponding bands show a very good 
overlap when qt ~ Mh- We thus conclude that, within the NLO theoretical uncertainty, the two 
calculations are perfectly compatible at the quantitative level in the large-gy region, qt ~ Mh- 

In Fig.|Sl(Fig.|ni) we plot the NLL+LO (NNLL+NLO) spectra for different choices of the resum- 
mation scale Q- We remind the reader that the resummation scale Q has to be chosen of the order 
of Mh- Variations of the resummation scale around Mh can be studied to estimate the uncertainty 
of the resummed calculation arising from not yet computed terms at higher logarithmic accuracy. 
In our quantitative study we consider four different values of Q, Q = 2Mh, Mh, Mh/2, Mh/^- 

We first comment on the behaviour at large transverse momenta, which is best visible looking 
at the plots on the right of Figs. |H1 and 13 We see that the NLL+LO cross section can become 
negative ii Q = 2Mh- This behaviour should not be regarded as particularly worrisome: it takes 
place when qt > Mh, where the use of the resummation formalism is not anymore justified. In 
general, the cross section has a better behaviour at large qt when the resummation scale has the 
values Q = Mh, Mh/^,, Mh/4:- In particular, at large-g^ the results of the fixed-order calculation 
at LO (NLO) accuracy are very well approximated by the NLL+LO (NNLL+NLO) calculation 
with Q = Mh/2; the line corresponding to the LO (NLO) results is not shown in the plot on the 
right of Fig. |H1 (Fig. since it is hardly distinguishable from the dotted and dot-dashed lines. 
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Figure 7: Comparison of the NLL+LO (left) and NNLL+NLO (right) hands with the NLO hand 
computed hy using NLO parton densities and 2-loop a^. 

The fact that the fixed-order behaviour at large qt is approximated better when Q is smaller is 
not unexpected. By varying Q, we smoothly set the transverse- momentum scale below which the 
resummed logarithmic terms are mostly effective; when Q is smaller, the resummation effects are 
confined to a range of smaller values of qt- 

To quantify the resummation-scale uncertainty on the cross section at small and intermediate 
values of qx-, we proceed as in the case of the renormalization and factorization scales, and we 
vary Q by a factor of 2 up and down from a reference value. We choose the reference value 
Q = because of the better quality of the behaviour of the corresponding cross section at 

large qt- From Fig. we see that at NLL+LO accuracy a scale variation between 1/AMh and 
Mh produces a variation of the cross section of about ±15% in the region around the peak. At 
NNLL-I-NLO accuracy (Fig.E)) the resummation-scale dependence is much reduced: when Q varies 
between Mh/4: and Mh the change in the cross section at the peak is about ±5%, i.e., smaller 
than the corresponding uncertainty from variations of the renormalization and factorization scales 
(see Fig. El). 

Throughout this section we used the MRST2004 set [70] of parton distribution functions at 
NLO and NNLO. The NLO and NNLO parton densities from Alekhin are currently being up- 
dated [7T]. The CTEQ [72] and GRV [72] groups do not include sets of NNLO parton densities. 
The parton distribution sets of MRST, Alekhin and CTEQ include estimates of experimental 
uncertainties, which lead to effects below to about 5% on the total cross section for Higgs boson 
production at the LHC. We do not expect significantly different results in the case of the qt 
cross section at the LHC, and we refer to Ref. ^7] for results and discussions about the effects of 
available parton densities on the total cross section. 
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Figure 8: NLL+LO spectra for different choices of the resummation scale Q at fixed jiji = 

The numerical results presented so far refer to the value Mh = 125 GeV of the Higgs boson 
mass. By varying Mh-, the typical features of the results are unchanged, the main difference 
being the decrease of the cross section as Mh increases. In Fig. we plot the NNLL+NLO 
spectra, normalized to the total cross section, for different values of the Higgs boson mass, Mh = 
125, 165, 200 and 300 GeV. For reference, the corresponding values of the NNLO total cross sections 
are (Tnnlo = 38.43,24.37,17.78 and 10.03 pb. As expected, the qx distribution becomes harder 
as Mh increases. The average value, (q't), of the transverse momentum increases almost linearly 
with increasing Mh-, and it is very roughly approximated by an effective lowest-order expression, 
(gr) ~Ca«s(M|)Mjj. 

The quantitative predictions presented up to now are obtained in a purely perturbative frame- 
work. It is known (see e.g. Ref. [IH] and references therein) that the transverse-momentum 
distribution is affected by non-perturbative (NP) effects, which become important as qx becomes 
small. In impact parameter space, these effects are associated to the large-6 region. In our per- 
turbative study the integral over the impact parameter turns out to be dominated by the region 
where 6 < 0.1-0.2 GeV~^, larger values of h being strongly suppressed by the resummation of the 
logarithmic terms in the gluon form factor. Thus we do not expect particularly-large NP effects 
in the case of Higgs boson production at the LHC. This expectation is in agreement with the 
findings in Refs. HOHlll . 

A customary way of modelling NP effects in the case of DY lepton-pair production is to intro- 
duce an NP transverse-momentum smearing of the distribution. This is implemented by multiply- 
ing the 6-space perturbative form factor by an NP form factor. Several different parametrizations 
of the NP form factor are available in the literature [UHllTlHTTj : the corresponding NP parameters 
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Figure 9: NNLL+NLO spectra for different choices of the resummation scale Q at fixed jiR = 

are obtained from global fits to DY data. 

In the case of Higgs boson production, the estimate of NP effects is obviously more uncertain, 
since we cannot exploit available experimental data. In Ref. [7H] we studied the impact of NP 
contributions on the qt spectrum of the Higgs boson, by applying the DY NP corrections of 
Refs. [TIHTH] to our resummed results at NLL accuracy. We also considered the effect of rescaling 
the DY NP coefficients by the factor Ca/Cf, to take into account the different colour charges of the 
initial-state partons {qq in the DY process, gg in Higgs boson production) in the hard-scattering 
subprocess. Alternatively, we used the NP coefficients extracted in Ref. 021 from a fit of data on 
T production, a production process that is more sensitive to the gluon content of the colliding 
hadrons. All these different quantitative implementations of NP corrections, although certainly 
not fully justified, can give an idea of the size of the NP effects on the Higgs boson spectrum. 

The results of Ref. (TH] show that the impact of the NP effects on the NLL resummed dis- 
tribution is definitely below 10% for g-r > 10 GeV, and it decreases very rapidly as qr increases. 
Moreover, when < 10 GeV, different parametrizations of the NP terms can lead to sizeably 
different relative effects, as a consequence of our present ignorance on the absolute value of the 
NP contributions. 

In view of these results, in the present paper we limit ourselves to considering a simple 
parametrization of the NP contributions. We multiply the 6-space resummed component W^(6, M, s) 
on the right-hand side of Eq. ()10|) by a NP factor, Snp, which includes a gaussian smearing of the 
form 

Snp = exp{-gNp . (89) 
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Figure 10: NNLL+NLO spectra for different values of the Higgs boson mass. The scales are set 
at the default value fJ^r = f^R = Q = Mh- 



The NP coefficient Qnp is varied in the range suggested by the study of Ref. Qnp = 1-67- 
5.64 GeV^. Note that this procedure, with these values oi gNP, well approximates the quantitative 
spread of NP effects found in Ref. [TH] at NLL accuracy. In Fig. ^2 we plot the effect of the NP 
smearing on our best perturbative predictions, as given by the results at NNLL+NLO accuracy. 
The inner plot shows the relative deviation from the NNLL+NLO perturbative result, as defined 
by the ratio 

_ da^^LL+NLO - daNNLL+NLO 

-r 5 lyuj 

O-^y NNLL+NLO 

where 'icrj^^^^_,_^^Q is the NNLL+NLO cross section, daNNLL+NLO, supplemented with the NP 
form factor. We see that the NP effects give deviations from the purely perturbative result that 
are below 10% for ^ 5 GeV. Comparing the inset plots in Figs. El and [TTl we also notice that 
the inclusion of higher-order contributions (going from NLL+LO to NNLL+NLO) and of NP 
contributions have a qualitatively similar effect at intermediate and small values of transverse 
momenta: both contributions make the distribution harder. At the quantitative level, A is much 
smaller than K — 1 when > 10 GeV, while A and K — 1 are comparable when qrp < 10 GeV. This 
points towards a non-trivial interplay between higher-order perturbative effects and NP effects at 
fixed value of the Higgs boson mass. 

In summary, the comparison of the NLL+LO and NNLL+NLO results from small (around 
the peak region) to intermediate (say, roughly, qj- <Mh/3) values of transverse momenta shows a 
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Figure 11: The NNLL+NLO perturbative results supplemented with the NP form factor in Eq. hSf^) . 
The upper (lower) curve at small qr is obtained with g^p = 1-67 GeV^ {dNP = 5.64 GeV^). 

nice convergence of the resummed QCD predictions for the qx spectrum of the Higgs boson at the 
LHC From this comparison and from the effects of variations of the renormahzation, factorization 
and resummation scales, we conclude that the perturbative QCD uncertainty of the NNLL+NLO 
results is uniformly of about 10% over this range of transverse momenta. The perturbative and 
NP uncertainty increases at smaller values of qx (see Figs. IHl and ITT|) : the perturbative uncertainty 
increases also at larger values of |^nM22j . The perturbative uncertainty on the NNLO cross 
section ^3], as estimated in the same manner (i.e. by comparing the NLO and NNLO results, 
and performing scale variations), is about 15% Our results on the qr spectrum are thus fully 
consistent with those on the total cross section, since the bulk of the events is concentrated at 
small and intermediate values of the Higgs boson qx- 



4 Conclusions 



In this paper we have considered the transverse-momentum spectrum of generic systems of high- 
mass M produced in hadron-hadron collisions. Following our previous work on the subject ISSIP, 
we have illustrated and discussed in detail a perturbative QCD formalism that allows us to resum 
the large logarithmic contributions in the small-g-r region {qx ^ M) and to consistently match 
the ensuing result to the fixed-order contributions in the large-g^ region {qx ~ M). The main 
features of our approach, that make it different from other implementations of 6-space resummation 
presented in the literature, are summarized below. 
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• The resummation is performed at the level of the partonic cross section. The parton dis- 
tributions are thus evaluated at the factorization scale fip, which has to be chosen of the 
order of the hard scale M. The resummation formula is then organized in a form that is in 
close analogy with the case of event shapes variables in hard-scattering processes j54H57j and 
threshold resummation in hadronic collisions [HHl [HH] : the various classes of logarithmic con- 
tributions are controlled by the QCD coupling as{fi%) evaluated at the renormalization scale 
fiR. This procedure naturally allows us to perform a systematic study of renormalization- 
and factorization-scale dependence, as is customarily done in fixed-order calculations. This 
should be compared with the other implementations of 6-space resummation, where the scale 
at which the parton distributions are evaluated is of the order of 1/6, which also necessarily 
requires an extrapolation of the parton distributions in the NP region. 

• The large logarithmic contributions are exponentiated in the form factor exp{^7v}, where 
the function Qj^f (see Eq. (HH)) is universal: it does not depend on the produced high-mass 
system, and it only depends on the flavour of the partons involved in the hard-scattering 
subprocess. More precisely (see Appendix various process-independent form factors 
control the various partonic channels. The process dependence, as well as the factorization- 
scale and factorization-scheme dependence is fully included in the hard-scattering coefficient 
Hn (see Eq. 

• We impose a constraint of perturbative unitarity through the replacement in Eq. ()16|): the 
6-space form factor exp{^jv(L)} is equal to unity at 6 = 0. This constraint has a twofold 
purpose. On one hand, it avoids the introduction of unjustified higher-order contributions 
in the small-6 region, which are present [ZH] in standard implementations of 6-space resum- 
mation. On the other hand, it allows us to recover the total cross section at the nominal 
fixed-order accuracy upon integration over qt- Note that, as a consequence, perturbative 
uncertainties at intermediate values of qt are reduced. 

The resummation formalism has been applied to the production of the SM Higgs boson in pp 
collisions. We combined the most advanced perturbative information that is available at present 
for this process: NNLL resummation at small qt and fixed-order perturbation theory at NLO 
at large qx- We developed a numerical code, named HqT that performs the calculation at 
NLL+LO and NNLL+NLO accuracy. In Sect. l3.ll we have presented a selection of results that can 
be obtained by our program at LHC energies. Owing to the unitarity constraint, the integral of our 
spectra at NLL-I-LO (NNLL+NLO) correctly reproduces the total NLO (NNLO) cross sections. 
The results show a high stability with respect to scale variations and an increasing stability when 
going from NLL+LO to NNLL+NLO accuracy. As summarized at the end of Sect. 13. this 
suggests that the uncertainty from missing higher-order perturbative contributions is under good 
control. 



A Appendix: Exponentiation in the multiflavour case 

In Sects. Em and 1231 we have discussed the exponentiation structure of the resummed component 
of the qt distribution. To simplify the notation and the presentation, we have limited ourselves 
to illustrating the case in which the partonic scattering involves a single fiavour of partons. This 
appendix is devoted to generalize the exponentiation to the case with partons of different fiavours. 
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To obtain Eq. (j43|) . the multiflavour analogue of Eq. (|12p. we start from the representation 
in Eq. (jl^ of the resummed partonic cross section W^^^ ^, and then we proceed as in Sect. 12.31 
The main difference with respect to the steps in Eqs. fl44|) - (jmj) is that the solution of the QCD 
evolution equations (Pn|l has the customary form"!" 

L/^(6>^Q2)=Pexp|y^^ , (91) 

where the symbol P on the right-hand side denotes the path ordering expansion of the exponential 
matrix. Because of its matrix structure, the exponential in Eq. ()91|) has only a formal meaning. 
To recast Eq. (j9T| in a true exponential form, we can perform a systematic logarithmic expansion 
of the solution of the Altarelli-Parisi equations, by using a well-known method that dates back, 
at least, to Ref. lED]. 

The evolution operator in Eq. ()91|) can be written in the following form [HU] (see also Ref. |HH 
for technical details): 

UM{hl/h\Q') = VM{as{hl/h'')) U^i;''\as{bl/b'),asm) ^^'(«s(Q')) , (92) 
where t/^^'' is determined by the lowest-order anomalous dimensions 7^^'', 

and the operator fulfils the following differential equation: 

^^f^ = 7^ 7^(«s) V^ias) + V^ias) ^ 7!^^ • (94) 
dlnas P(as) Po 

The evolution equation (jH!^ can be solved by diagonalizing the anomalous dimensions matrix 
7^'*, which has three different eigenvalues Ji^lj'- one eigenvalue in the flavour non-singlet sector 
{i = NS), and two eigenvalues in the flavour singlet sector {i = ±). The solution of Eq. is 



C/r'(as(6o/''').as(g')) 



=NS,± 



' E^S , (95) 



WsiPl/h' 



where E^j denotes the projector onto the flavour eigenspace corresponding to the eigenvalue 7J■^]^r. 
By inspection of Eq. (jMj) . we see that it can be solved by performing a perturbative expansion. 



00 



TT 
n=l 



and the perturbative coefficients vj^^ are obtained in terms of the anomalous dimensions coeffi- 
cients 77v^^'* and the (3 function coefficients jSk with k < n. For example, the first-order coefficient 
V^^^ is given by 



Tin this appendix we use the boldface notation X to denote the flavour space matrix whose matrix elements 
are Xab = (X)^^. 
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We now come back to the right-hand side of Eq. (P^ . The evolution operator UN{bl/b'^, fj,p) 
is rewritten as U Nibl/b^ , fi^) = U Nipl/b'^ .Q"^) Uj^iQ'^, /J-f)- Then U Nipl/b'^ ,Q'^) is replaced by 
the expression in Eq. (jU^ . Equation (ji^ thus becomes 

y\^^,n{b,M■a^{^^l),^^l,^^l) = J2^T,fMM'),M) H[{as{M')) 

c 

5^ [y^l(as(Q^))C/^(g^/.^)]^^^ [V],\asm)U^{Q',l^l 



X 



62 b 



12,62 



X 



ai,bi 

where we have defined the perturbative function 



(9^ 



(99) 



71=1 



and inside the curly brackets we have collected all the factors, Sc, C ^ and t/j^^'' , that depend on 
the impact parameter b. These factors contain the logarithmically-enhanced contributions that 
have to be resummed and organized in exponential form. The factor Sc can be rewritten as 



S,{M,b) = S,{M,bo/Q) exp{g,{as{fil),L-Myfi%MyQ^)} 

where (see Eq. 



g,{as{fi'n),L;Myfi%MyQ'' 



The factor C/5J°^ is 



dq 



2 r 



A,(as(g')) In— + i?,(as(g')) 



(100) 



(101) 



j=NS,± 



(102) 



where (see Eq. 

g,,^{as{fil),L;Myf,l,MyQ') 



li,N 



Po as{bl/b'^) I3q Jqn,2 q 



K /?(«s(g^)) 



(103) 



The factor Cca,N can be written as 

Cca,N{as{bl/b^)) = C,a,N{as{Q^)) exp{g,a,N{as{fil),L;Myfi%MyQ^)} 
where (see Eq. (jSZj)) 



gca,N{as{fiji),L;Myfi%MyQ^' 



/?(«s(g')) 



d\n.Cca,N{as{q^)) 
(ilnas(g2) 



(104) 



(105) 
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Note a key point: Eq. (jlU4j) does not regard the flavour matrix Gat, but rather its matrix element 
Cca,N- Therefore, its right-hand side involves a true c-number exponential instead of a formal 
matrix exponential. 

Inserting Eqs. (jlOOj) . p02j) and p04|) in Eq. ()98|) . we eventually obtain the final exponentiated 
result in Eq. (jlBI), namely 

{1} 

X exp{e?|,},;v(«s(/i?j),i^;MV/i^,MVQ2)| ^ 

where the sum extends over the following set of flavour indices: 

{/} = c,c,ij,ai,bi . (106) 

The exponent Q{i}^n of the universal form factor and the pro cess- dependent hard factor 'H\^ 
are 

'3{I},N = Qc + '3i,N + Qj,N + Qcai,N + Gchi,N , (107) 

X S,{M,bo/Q) C,a„N{as{Q')) 5a,^(as(Q')) (108) 



X 



J bib 



From Eqs. pUl|) . ()103p and pOSp we see that Q{i},n in Eq- P07|) has exactly the integral repre- 
sentation of Eq. (fT9|l . The logarithmic expansion (see Eq. (fT4|l ) of and ^j^A? starts at LL and 
NLL accuracy, respectively. The term Qca,N starts at NLL accuracy in the flavour off-diagonal 
case (c 7^ a) and at NNLL accuracy in the flavour diagonal case (c = a). The hard function 
Ti-ab does not depend on the impact parameter b. It can be perturbatively expanded in powers 
of as(/^|;) (with fiji ~ M), since the various factors on the right-hand side of Eq. pOSj) involve 
only scales {M,Q, fip) that are of the order of the hard-scattering scale M. 

We conclude this appendix with a comment on the solution (j^^ of the Altarelli-Parisi evolution 
equations and its relation with the resummation in Eq. (fT^ of the logarithmic contributions to 
the impact-parameter form factor exp{Q^iy iy{as, L)}. 

The evolution operator UN{bl/b'^,Q'^) does not contribute to the LL function g^^\asL) in 
Eq. (jl4j) . It starts to contribute to the resummation at the level of the NLL function g^^\a^L). 

Indeed, from Eqs. (jl(J2p and (jl(J3p we see that t/^°''(as(^o/^^)' Q^s(Q^))) the solution of the evolu- 
tion equations at the lowest-perturbative order, contributes to the NLL terms a^L^. The higher- 
order corrections to the evolution equations are taken into account by the operator VAr(as) in 
Eq. ()92|) . The role of these corrections can be examined by organizing them in classes of logarith- 
mic contributions aKas-Z^)". Using the expansion in Eq. (jHSl), Eq. (jU^ gives 

c/^(&^/&',Q')=c/rH«s(fe^/&'),«s(g'))+0(«r'^"+') , (^>o), (109) 

c/iv(&^/&^Q') = c/^°H«s(&^/&'),«s(Q')) + (^(«r'^"+') . (^>o), (110) 
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where U^^^'^^ is the customary solution (HHll^ of the evolution equations at NLO: 



)) 



"^(^')t/rK(&^/&^),«s(Q^))< . (Ill) 



TT 



The terms denoted by C(ag+^L"+^) on the right-hand side of Eq. (fTnn|l contribute at NNLL 
accuracy (they are of the same logarithmic accuracy as those in the function a^g'^^^a^L) in 
Eq. (HH)). Analogously, the terms denoted by C(ag"^^L"+^) on the right-hand side of Eq. ()109|) 
contribute at NNNLL accuracy (they are of the same logarithmic accuracy as those in the function 
alg^ff{asL) in Eq. (HH)). Therefore, to resum the NLL (NNLL) contributions to the form factor is 
sufficient to implement the solution of the evolution equations at the LO (NLO). Note, however, 
that, to be consistent with the resummed logarithmic expansion, the scale dependence of the 
running couplings as(^o/^^) '^s(Q^) in Eq. p09|) (Eq. piO|) ) has to be evaluated at the NLO 
(NNLO). 



B Appendix: Bessel transformation of logarithmic contri- 
butions 

This appendix is devoted to the computation of the Bessel transformation of logarithmic contri- 
butions. 

We recall the definition of the the functions In^qr/Q) introduced in Eq. (ffljl : 

U^t/Q) = db - Mbqr) In" + 1 j . (112) 

These integrals are easily evaluated in terms of derivatives of the corresponding generating function 
7(x;e): 

d 



J„(x) = Ihn (^-J J(x;e), (113) 
where 

00 ^ 

n=0 

Inserting Eq. (jll2p in the right-hand side of Eq. (jll4j) , we have 

I{x-t) = dt-Joitx) + , (115) 

and this integral can be expressed jH21 as follows in terms of K^{x), the modified Bessel function 
of imaginary argument (see Eq. (fTSj) ): 
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Inserting Eq. in Eq. (jll3|) and using the relation 

CXD ^ 

r(l-e)=exp{7Ee + 5^^ae'=} 



(117) 



k=2 



where Cn is the Riemann zeta-function (C2 = 7r^/6 = 1.645 . . . , ^3 = 1.202 . . . ), the integrals 
can straightforwardly be expressed in terms of the derivatives, K'f'\z), of the Bessel function with 
respect to its index z/: 



These derivatives have the following integral representation: 



poo 

Jo 



2n ^—zcosht 



cosht 



{2n+l)^^N, _ 2n + 1 /"°° ^2n g-2cosht 



(118) 

(119) 
(120) 



which can simply be obtained from Eq. (j75p . 



As discussed in Sect. 12.41 the computation of the finite component of the qt distribution 
requires the evaluation of the functions /„(a;) when x > 0^. In particular, the computation up to 
NLO (see Eq. (fTHjl) requires /„(x) with n = 1,2,3,4; these functions are 



X 



hix) 



2bo 



X 



Ki{box) Inx - K[^\box) 



7 ( \ 

h[x) = 

X 



Ki{box) (In^x - C2) - 2K[^\box) \nx + K^ibox) 



'(2), 



(121) 

(122) 
(123) 



h{x) = — 

X 



Ki (box) (In^ X - 3C2 In X + 2(3) - 3K[^^ (box) (In^ x - C2) 



+ 3K[^\box) Inx- Kl-^'ibox) 



(3), 



(124) 



The functions /„(x) diverge when x 0. To examine the divergent behaviour at small values 
of X, we introduce the functions In{x) and the corresponding generating function I{x; e): 



In{x) = In{x) [1 + 0{X' 



(125) 



I{x;e) = l{x;e)[l + 0{x^)] . (126) 

Using the small-x behaviour of the Bessel function Ki+^{x) [H21 and performing the small-x limit 
of Eq. jnH), we get 

7{x;e) = -eD{e) (x^)"^"^ , (127) 



■I- The behaviour of In{x) when x = is discussed at the end of this appendix. 



41 



where 



D(e) 



2 V' Til 



r(l-e) 



exp 



k=l 



2k + 1 



J2k+l 



(128) 



Note that the functions In{x) exactly correspond to the following Bessel transformations: 



Iniqr/Q) = 



b fO%'^ 
db-MbqT) In" PV 
^ V "0 



(129) 



as can be checked by performing the limit gj- ^ of Eq. (jll2|) or by verifying that the generating 
function in Eq. (jl27p has the following integral representation: 



I{x; e) = - 6( 



-2e 



dt t^^^' Jo{tx) . 



(130) 



The relation between In^qx/Q) and the small-gr limit of Inigr/Q) is not unexpected in view of the 
discussion in Sect. 12.21 The integral in Eq. (jll2p originates from Eq. (jl29p after the replacement 
L = ln{Q'^b'^ /bl) — > L = ln(l + Q'^b'^/b'^) at the integrand level: when qt — > 0, such a replacement 
has no effects on the singular behaviour at any logarithmic accuracy. 

Though In{.<iT/Q) and Inilr/Q) coincide when — 0, they behave quite differently at very 
large values of qr- When a; ^ 00, from Eqs. ()116|1 and ()127|1 we get 



In{x) 
7„(x) 



in— 1 



n 



In 



box ^ 

2 

X 1 + 



1 + 



1 



Inx 



n-l 



X^ 



\nx 



(131) 
(132) 



Note, in particular, that /„(a;) is integrable over when x —>■ 00, whereas /n(x) it is not. 



The function J„(x) can easily be computed by performing the n-th derivative of the generating 
function p27j) with respect to the parameter e. To present the result, we first exclude the singular 
point X = and consider only the region x > 0. Since the generating function depends on x only 
through the factor (x^)"^"*^, x^ In{x) is simply a polynomial of degree n — 1 in the variable Inx^: 



/„(x) = - 



^ n— 1 



x^ — ' k\{n — k — 

fc=0 ^ ' 



4 In"-'^-^ 4 



X 



2 ' 



X > 



(133) 



where the coefficients dn are obtained from Eq. 



dr. 



de 



Die] 



J e=0 



The value of the first few coefficients is 

do = 1, di = d2 = 0, 4 = -4(3, 



0, 4 = -48(5, de 



160CI 



(134) 



(135) 
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The result in Eq. (|133|) agrees with that in Ref. jHSl, where one can find the numerical values of 
dn with n < 19 [dn = 2"'bn{oo), where bn{oo) are given in Table 1 of Ref. [12] )• The small-x limit 
of Eqs. (IT^ - (IT^ thus gives 



x"^ 



(136) 

Note that, since c^i = = (^4 = 0, In{x) can be expressed in a simple form to very high logarithmic 
accuracy. For example, we have: 

n f- n-i 1 2^ {n-l)\ 1 2^ {n - l)\ ^^„_6 ^ 



1 x^ 3 (ra — 4)! x^ 5 (n — 6)! x^ 



2 .2 (n-1)! 7 1 , o f 1 



X > . 



(137) 



We now discuss how to deal with the region around the singular point x = 0. We first split 
the X range in a large-x (x > Xq) and a small-x (x < Xq) region, where the parameter xq can be 
chosen arbitrarily. Setting xq = 1, we have 



J„(x) = J„(x) e(x - 1) + /„(x) 6(1 - x) . 



'138) 



In the large-x region, which excludes the point x = 0, In{x) is given by Eq. ()133|) . In the small-x 
region, to properly treat the singularity at x = 0, we have to consider the generating function in 
Eq. (|127|) and use the expansion 



(x^)-^-^ 0(1 - x) = -- 5(x^) + 



X 



2\-e 



X^ 



6{x') + y: 



n=0 



X^ 



x^ 



(139) 



where the plus-distribution is customarily defined by its action onto any function hlx"^) that is 
finite at X = 0: 



dx'^ /i(x^ 



x^ 



x^ 



x^ 



x^ 



(140) 



Therefore the generalization of Eq. p33j) to include the point x = is 

rt-l 



/„(x) = dn (5(X^) - ^ 



k=0 



k\{n-k- 1)! 



di 



x^ 



x^ 



< X < 1 . 



(141) 



The procedure described in Eqs. ()138|1 and ()139|1 can also be applied to properly define the 
integrals /n(x) around the point x = in the small-x region. Choosing xq = oo, the final result is 
equivalent to start from /„(x > 0), the expression of /„(x) when x 7^ (for example, Eqs. p21|l - 
fjl24|) ). and then introduce a generalized plus-prescription that acts in the entire range < x < 00. 
Formally we can write 



IJx) 



In{x > 0) 



where the generalized plus-distribution is defined as 



dx^ h{x^) /„(x > 0) 



dx"^ [h{x^) - h{0)] In{x > 0) . 



(142) 



(143) 
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The choice Zq = oo to define the plus-prescription in the case of J„ is feasible since In{x) (unlike 
In{x)) is integrable over when x ^ oo. This choice simplifies the definition of J„ since the 
right-hand side of Eq. ()142|) (unlike Eq. ()141|) ) does not contain any contact term proportional to 
S{x'^). The contact term vanishes since the integrand factor ln"(l + Q'^h'^ /h"^) in Eq. ()112|) vanishes 
at 6 = 0. The vanishing of the contact term is thus ultimately related to the unitarity constraint 
in Eqs. (|HD and (HH). 
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